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^trailing  edge  for  the  flow.  With  this  model,  the  slip  line  boundaries  create 
an  exponential  damped  sinusoidal  relaxation  of  the  disturbances  as  compared 
to  algebraic  decay  for  analogous  unconfined  flows.  A  similar  observation 
applies  to  wind  tunnel  far  fields.  Ihe  nonlinear  case  differs  from  the  lin¬ 
ear  one  in  that  the  amplitude  of  the  downstream  propagated  disturbance  in¬ 
teracts  nonli nearly  with  the  near  field  in  the  former  case.  For  these  wall 
jets,  acceleration  to  criticality  is  accomplished  by  stream  tube  contrac¬ 
tions  and  throats  induced  by  upstream  influence  of  the  turning.  In  the  co- 
flowing»case,  selection  rules  defining  various  jet  flow  regimes  in  terms  of 
the*jet  transonic  similarity  parameter  and  pressure  level  are  given.  Condi¬ 
tions  are  prescribed  which  define  the  penetration  length  for  supersonic  dis¬ 
turbances  downstream  from  the  exit  as  well  as  the  location  of  the  throats  in 
the  jet.  The  criticality  of  the  flow  far  downstream  is  also  quantified  in 
terms  of  these  parameters.^-. 

Models  are  described  for  tlT5~ complete  problem  for  a  tangentially  blown 
a’rfoil  in  incompressible  and  transonic  flow.  An  important  feature  of  these 
configurations  is  the  foregoing  wall  jet  structures.  A  systematic  theory 
based  on  matched  asymptotic  expansions  for  the  fine  structure  of  the  jets 
away  from  their  exits  is  developed  which  shows  that  in  a  small  deflection 
approximation  the  pressure  jumps  associated  with  Spence's  theory  prevail, 
even  if  the  flow  is  rotational  and  compressible.  Furthermore,  the  asympto¬ 
tic  developments  provided  allow  further  systematic  refinements  in  which 
the  effects  of  initial  skewness  and  vorticity  inaccessible  to  other  theories 
can  be  assessed.  A  study  of  the  trailing  edge  region  of  the  blown  airfoil 
problem  reveals  that  for  finite  trailing  edge  angles,  the  dividing  stream¬ 
line  leaves  tangent  to  the  higher  stagnation  side  if  upper  and  lower  stag¬ 
nation  pressure  are  unequal.  If  these  pressures  are  equal,  the  dividing 
streamline  leaves  along  the  bisector.  Computational  results  based  on  an 
inviscid  model  reveal  that  significant  enhancements  of  lifting  pressures 
are  possible  with  tangential  blowing  on  transonic  airfoils.  Comparisons 
with  experiment  indicate  the  need  for* refinements  incorporating  wave  inter¬ 
action  phenomena  near  the  jet  exit  as  well  as  viscous  interaction  processes 
in  the  downstream  portion  of  the  wall  jet.  Based  on  this  requirement,  a 
finite  difference  module  has  been  developed  which  characterizes  the  develop¬ 
ment  of  blown  laminar  boundary  layers  and  wakes  in  the  transonic  reg-ime. 
Solutions  arising  from  this  phase  are  being  utilized  to  quantify  not  only 
viscous  interaction  effects  over  blown  airfoils,  but  also  the  degree  of 
delay  in  shock-induced  separation  achievable  with  slot  injection.  Sub¬ 
stantial  downstream  movement  in  the  shock-induced  separation  point  has 
been  inexpensively  quantified  with  the  computational  procedure  as  compared 
to  experimental  methods.  The  viscous  studies  show  also  that  the  propulsive 
wake  generated  by  blowing  gives  a  nearly  Gaussian  normalized  velocity 
profile,  in  spite  of  a  lack  of  self  similarity  which  occurs  in  submerged 
and  self-propelled  cases.  This  result  is  of  significance  to  assessment 
of  net  thrust  in  experimental  simulations  and  theoretical  far  field 
modeling. 
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SUMMARY 

A  study  of  transonic  wall  jets  and  tangentially  blown  wings  has  been 
performed  by  using  asymptotic  and  computational  methods.  For  the  portion  of  the 
effort  dealing  with  wall  jets,  both  submerged  and  coflowing  cases  have  been 
Investigated.  For  the  submerged  configurations,  the  nature  of  the  decay  process 
for  disturbances  emanating  from  the  jet  exit  has  been  examined.  In  this 
analysis,  the  Kutta  condition  on  the  nozzle  rim  can  be  satisfied  merely  by 
requiring  streamwise  continuity  of  the  potential  across  the  rim  as  a  trailing 
edge  for  the  flow.  With  this  model,  the  slip  line  boundaries  create  an 
exponential  damped  sinusoidal  relaxation  of  the  disturbances  as  compared  to 
algebraic  decay  for  analogous  unconfined  flows.  A  similar  observation  applies 
to  wind  tunnel  far  fields.  The  nonlinear  case  differs  from  the  linear  one  In 
that  the  amplitude  of  the  downstream  propagated  disturbance  Interacts 
nonllnearly  with  the  near  field  In  the  former  case.  For  these  wall  jets, 
acceleration  to  criticality  Is  accomplished  by  stream  tube  contractions  and 
throats  induced  by  upstream  Influence  of  the  turning.  In  the  coflowing  case, 
selection  rules  defining  various  jet  flow  regimes  in  terms  of  the  jet  transonic 
similarity  parameter  and  pressure  level  are  given.  Conditions  are  prescribed 
which  define  the  penetration  length  for  supersonic  disturbances  downstream  from 
the  exit  as  well  as  the  location  of  the  throats  in  the  jet.  The  criticality  of 
the  flow  far  downstream  Is  also  quantified  in  terms  of  these  parameters. 

Models  are  described  for  the  complete  problem  for  a  tangentially  blown 
airfoil  in  incompressible  and  transonic  flow.  An  important  feature  of  these 
configurations  is  the  foregoing  wall  jet  structures.  A  systematic  theory  based 
on  matched  asymptotic  expansions  for  the  fine  structure  of  the  jets  away  from 
their  exits  is  developed  which  shows  that,  in  a  small  deflection  approximation, 
the  pressure  jumps  associated  with  Spence's  theory  prevail,  even  if  the  flow  is 
rotational  and  compressible.  Furthermore,  the  asymptotic  developments  provided 
allow  further  systematic  refinements  in  which  the  effects  of  initial  skewness 
and  vorticity  inaccessible  to  other  theories  can  be  assessed.  A  study  of  the 
trailing  edge  region  of  the  blown  airfoil  problem  reveals  that  for  finite 
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trailing  edge  angles,  the  dividing  streamline  leaves  tangent  to  the  higher 
stagnation  side  if  upper  and  lower  stagnation  pressures  are  unequal.  If  these 
pressures  are  equal,  the  dividing  streamline  leaves  along  the  bisector. 
Computational  results  based  on  an  inviscid  model  reveal  that  significant 
enhancements  of  lifting  pressures  are  possible  with  tangential  blowing  on 
transonic  airfoils.  Comparisons  with  experiment  Indicate  the  need  for 
refinements  incorporating  wave  interaction  phenomena  near  the  jet  exit  as  well 
as  viscous  interaction  processes  in  the  downstream  portion  of  the  wall  jet. 
Based  on  this  requirement,  a  finite  difference  module  has  been  developed  which 
characterizes  the  development  of  blown  laminar  boundary  layers  and  wakes  in  the 
transonic  regime.  Solutions  arising  from  this  phase  are  being  utilized  to 
quantify  not  only  viscous  interaction  effects  over  blown  airfoils,  but  also  the 
degree  of  delay  in  shock  induced  separation  achievable  with  slot  injection. 
Substantial  downstream  movement  in  the  shock-induced  separation  point  has  been 
inexpensively  quantified  with  the  computational  procedure  as  compared  to 
experimental  methods.  The  viscous  studies  show  also  that  the  propulsive  wake 
generated  by  blowing  gives  a  nearly  Gaussian  normalized  velocity  profile,  in 
spite  of  a  lack  of  self  similarity  which  occurs  in  submerged  and  self-propelled 
cases.  This  result  is  of  significance  for  assessment  of  net  thrust  in 
experimental  simulations  and  theoretical  far  field  modeling. 
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NOMENCLATURE 

a  Characteristic  submerged  jet  height 

b  Slip  line  deflection,  b  coefficient  in  boundary  layer  momentum 
equation 

c  Density  ratio  pe/p  in  boundary  layer,  chord  length 
d  Jet  exit  height  In  physical  plane 

f  wall  and  airfoil  shape  function,  similarity  factor  for  stream 
function 

g  Jet  displacement  function 

h  Metric  coefficient 

t  Unit  vector  in  Cartesian  x  direction 

J  U.,it  vector  In  Cartesian  y  direction 

m  Ratio  of  kinetic  energy  per  unit  mass  to  freestream  stagnation 
enthal py 

m.j  1  =  Integer,  coefficients  In  boundary  layer  momentum  equation 
p  Perturbation  pressure,  Fourier  transform  variable 

q  Velocity  vector 

q  Speed,  dynamic  pressure 

r  Polar  radius 

Sj  Arc  length  along  airfoil  or  wall 
p  Velocity  component  In  x  direction 

\i'  Velocity  perturbation  about  freestream 

v  Velocity  component  In  y  direction 

w  Complex  velocity 

x  Cartesian  or  curvilinear  coordinate  in  freestream  direction  or 

along  parallel  to  wall 

x*  Cartesian  coordinate  in  freestream  direction 
y  Cartesian  or  curvilinear  coordinate  normal  to  freestream  or  to 
wall 

y*  Cartesian  coordinate  normal  to  freestream  direction 
z  Complex  variable 
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A  Peak  velocity  parameter,  constant  in  empirical  fit  for  external 
pressure  in  wake  region,  constant  in  nose  shape 
B  Constant  in  empirical  fit  of  external  pressure  in  wake  region, 
constant  in  nose  shape 
Cj  Blowing  coefficient 

C  Chapman  constant 

Cpp  Far  field  constant 

CL  Lift  coefficient 

Cp  Pressure  coefficient 

D  Reduced  exit  height 

E  Normalized  stagnation  enthalpy  ratio 

F  Wall  slope  function,  complex  potential 

G  Green's  function,  reduced  jet  velocity  function,  pressure  gradient 

parameter,  E' 

H  Heaviside  function,  stagnation  enthalpy 

J  Jet  momentum  efflux 

K  Transonic  small  disturbance  parameter 

L  Wall  length 

M  Mach  number,  moment  of  momentum 

Freestream  Mach  number 
N  Number  of  supersonic  points 

P  Pressure 

R  Radius  of  curvature,  reflection  coefficients 
Re  Reynolds  number 

S  Slip  line  implicit  function  for  shape 

T  Karman  Guderley  operator 

U  Freestream  velocity 

W  Typical  flow  quantity 

a  Angle  of  attack,  reflection  coefficient 

£  /m£-T,  Prandtl  Glauert  factor,  angle  between  dividing  streamline 

and  tangent  to  lower  surface 
y  Specific  heat  ratio 
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6  Characteristic  flow  deflection 

6*  Displacement  thickness 

t  Eddy  viscosity 

C  Boundary  layer  similarity  coordinate 

n  Dumnty  coordinate  In  y  direction.  Intermediate  variable, 

curvilinear  coordinate  normal  to  wall 
e  Momentum  thickness,  polar  angle 

\  Mach  number  ratio  appearing  in  reflection  coefficient 
u  Viscosity 

i  =  1,2,3  quantities  appearing  in  energy  equation 
v  Trailing  edge  angle 

5  Dummy  coordinate  In  x  direction,  curvilinear  coordinate  along 

parallels  to  wall 
p  Density 

t  Characteristic  slope  of  jet  free  boundaries 
$  Perturbation  velocity  potential  In  external  flow 

V  Stream  function 

ui  Coflowing  jet  overpressure  parameter,  angle  between  dividing 
streamline  and  upper  surface,  vorticity 
a  Laplacian 

r  Circulation 

$  Reduced  velocity  potential 


Subscripts 

e 

u 

i 

j 

FF 

SL 

00 

5 


Refers  to  edge  of  boundary  layer 
Upper 

Lower 

Jet 

Far  field 
Slot 

Freestream 

Component  in  5  direction 
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n  Component  In  n  direction 

0  Zeroth  order  quantity 

1  First  order  quantity 

Superscripts 

+  Above  dividing  streamline 

Below  dividing  streamline 

Special  Symbols 

<  >  Average 

[  ]  Jump 

~  Overbar,  refers  to  average  of  quantity  under  it. 
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1.0  INTRODUCTION 

Tangential  blowing  to  avoid  the  effects  of  separation  on  airfoils  and  to 
produce  superclrculatlon  lift  enhancement  has  received  extensive  consideration 
at  low  speeds.  A  natural  application  of  these  capabilities  to  current  tactical 
requirements  associated  with  runway  denial  and  ship  dispersal  is  of  great  mili¬ 
tary  Interest.  These  scenarios  strongly  suggest  the  value  of  V/STOL  capability. 
Predictive  models  for  the  aerodynamic  characteristics  of  jet-flapped  and  upper 
surface  blown  wings  associated  with  this  application  have  been  primarily  addressed 
through  linearized  surface  singularity  and  conformal  mapping  methods.1-4  Sys¬ 
tematic  and  mathematically  consistent  procedures  integrating  viscous  effects  in 
such  a  calculation  have  received  only  limited  attention,  and  appear  not  to  have 
addressed  the  question  of  the  necessary  blowing  requirements  to  avoid 
separation. 

At  transonic  speeds,  tangential  blowing  has  been  proposed  as  a  means  of 
mitigating  the  adverse  aerodynamic  consequences  of  shock  boundary  layer 
interaction.  In  particular,  the  quest  for  high  maneuverability  in  this  Mach 
number  range  suggests  enhanced  buffet-free  lift  through  supercirculation  and  jet 
momentum  effects.  Quantification  of  these  gains  as  well  as  others  associated 
with  such  non-aeronautlcal  applications  as  circulation  controlled  rotors, 
diffusers,  convective  film  cooling,  and  chemical  laser  aerodynamic  windows  requires 
further  knowledge  of  the  structure  of  the  associated  wall  jets  and  the  nonlinear 
interaction  with  the  surrounding  flow. 

Toward  obtaining  theoretical  models  which  could  describe  these  phenomena, 
ONR  has  had  a  research  program,  "Transonic  Wall  Jets  and  Upper  Surface  Blown 
Mings,"  underway  since  January  1976.  During  the  course  of  this 
investigation,  analytical  and  computational  procedures  have  been  utilized  to 
gain  further  insight  into  tangentially  blown  airfoil  flows.  The  work  has 
encompassed  detailed  investigations  of  the  wall  jets  and  their  relationship  to 
the  aerodynamic  augmentation  of  the  blown  airfoils.  By  contrast  to  the  surface 
singularity  procedures  utilized  in  the  other  early  studies  of  blown  wings,  the 
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effort  to  be  described  has  utilized  computational  methods  which  are  capable  of 
handling  nonlinear  phenomena  and  viscous  interactions  that  occur  in  the  transonic 
regime  in  a  mathematically  consistent  manner. 

This  report  summarizes  these  studies  and  provides  recommendations 
for  future  research.  It  is  based  on  a  number  of  papers  which  have  been 
disseminated^-?  -jn  connection  with  the  contract  and  other  resu1*s  which  will  be 
published  in  the  open  literature  in  the  future.  The  first  portion  of  this 
report  summarizes  the  early  work  of  the  contract  dealing  with  transonic  wall 
jets.  In  subsequent  sections,  a  study  of  inviscid  phenomena  over  incompressible 
and  transonic  blown  wings  is  provided.  Later  portions  describe  recent  work 
in  viscous  effects  dealing  with  delay  of  separation  and  the  structure 
of  wakes  downstream  of  propulsive  wings. 
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2.0  A  STUDY  OF  INVISCID  SUBMERGED  TRANSONIC  WALL  JETS 

An  essential  element  of  the  flow  field  over  propulsive  wings  such 
as  those  employing  upper  surface  blowing  is  transonic  wall  jets.  In  such 
phenomena,  a  transonic  primary  flow  or  wall  jet  entrains  an  ambient  or  coflowing 
stream  through  turbulent  mixing  processes. 

Existing  models  for  such  wall  jets  stress  the  incompressible  treatment 
of  these  phenomena  using  eddy  viscosity  and  energy  methods.  Correspondingly, 
there  is  a  need  for  simulations  that  include  the  effects  of  nonlinearities,  mixed 
flow,  and  wave  Interactions  on  the  development  of  the  wall  pressure 
distributions  and  overall  augmentation  forces. 

Previous  Investigations  of  related  phenomena  are  limited  to  the 
treatment  of  inviscid  shockless  free  jets,  and  include  the  work  of  Chaplygin,® 
Frankl  ,9  and  Guderley10,  all  of  which  employ  hodograph  methods.  To  study  shock 
development  and  mixed  flow  phenomena,  we  have  applied  modern  relaxation  methods 
to  treat  arbitrary  jet  exit  velocity  distributions  and  assess  the  influence  of 
an  adjacent  wall  boundary. 

In  this  section,  the  computational  model  will  be  discussed  from 
analytical  and  numerical  viewpoints  regarding  nonlinear  inviscid  phenomena.  In 
analogy  to  unbounded  cases  such  as  airfoil  flows,  the  far  field  is  employed  to 
condition  the  numerical  problem  and  provide  useful  information  about  the  decay 
of  disturbances.  Both  free  and  wall  jets  are  discussed  for  several  examples 
illustrating  various  features  of  this  class  of  flows. 

2.1  Formulation 

Referring  to  the  physical  configuration  depicted  in  Fig.  1,  a  jet  is 
shown  exhausting  from  the  exit  0C  bounded  by  the  wall  0Q  and  a  mixing  layer 
which  has  been  idealized  as  the  slip  line  CB.  This  approximation  neglects 
turbulent  diffusion  processes  in  the  study  of  wave  interactions  with  the  shear 
layer,  but  these  phenomena  can  be  incorporated  in  later  refinements. 

Furthermore,  it  will  be  assumed  that  wall  and  jet  turning  angles  are  small.  In 
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contrast  to  the  usual  jet  formulations,  ii  which  an  upstream  cowl  shape  is 
specified,  or  stagnation  conditions  are  assumed,8”12  this  analysis  will  treat  a 
specified  exit  Mach  number  distribution.  Additional  assumptions  are 
Irrotatlonallty  and  subsonic  conditions  Infinitely  far  downstream.  The  methods 
applied  here  can  be  generalized  to  cases  where  these  restrictions  are  not 
present.  Finite  length  walls  OQ  are  considered  In  keeping  with  relevance  to 
upper  surface  blown  wings  and  other  propulsive  lift  devices. 

Returning  to  Fig.  1,  the  equations  of  slip  lines  Sj(x,y),  S2(x,y),  and 
the  wall  boundary  B(x,y)  are  assumed  as 


CB:  Sj  =  y -  d  -  dGjfx)  =  0 
QA:  S2  =  y  -  6G2(x)  ~  G 
OQ:  B  *  y  -  $f(x)  =  0 

where  6  is  a  characteristic  flow  deflection  parameter.  In  a  small  disturbance 
limit,  in  which  the  scaled  jet  exit  height  D  =  ds1^3,  the  wall  length  L,  and  the 
transonic  similarity  parameter  K  =  (1  -  M2)/62^3  are  held  fixed,  6  0,  the 

oo 

asymptotic  expansions  of  the  velocity,  pressure  P,  and  density  p  are* 

I  (x.yjM^.d.d.L)  =  [1  +  <52/34>x(x,Y;K,D,L)  +  ...]t  +  [6*  +  ...If  (1.1a) 


P/P  =  1  -  Y«  '  *¥  +  ... 

00  X 


p/p  -  1  +  6 2^30  +  ... 


(1.1b) 


(l.lc) 


where  the  subscript  ~  signifies  conditions  at  x  =  «,  $  is  the  perturbation 

o 

potential,  P^  is  the  ambient  pressure,  U  =  a  M  ,  a  =  vP  /p  ,  p  is  the 
00  00  00  00  00  00  00 

density,  a^  is  the  speed  of  sound,  §  is  the  flow  velocity,  and  the  scaled 

1/3 

coordinate  Y  =  6  y  is  also  fixed  in  the  limit.  If  a  further  transformation  is 
introduced  in  which  y  =  /K  Y,  a  boundary  value  problem  can  be  formulated  for  the 
case  of  an  elliptic  far  field.  The  region  can  be  considered  as  the  rectangular 

*Equations  are  numbered  consecutively  from  the  beginning  of  each  section. 
References  to  equations  in  other  sections  provide  the  section  number  of 
the  referenced  equation. 
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domain  shown  in  Fig.  2,  corresponding  to  transfers  of  the  boundary  conditions  to 
the  appropriate  undisturbed  streamlines  allowed  by  the  small  disturbance 
limit.  Dropping  the  tildes,  the  following  small  disturbance  equation  holds 
inside  OQABCO 

A<|>  =  (a2/3x2  +  32/ay2 )  4,  =  (y+l)(3U2/ax)/2K  ,  (u  =  ^x)  (1.2) 

where  we  define  function  u  and  $  such  that 


u(x,y)  =  u,(x,y) 
<Ji(x,y)  =  <s(x,Y) 


Invoking  continuity  of 

with  a  =  D/K 

pressure  and  flow  tangency  along 

the  slip  lines,  we  have 

<|>(x,a)  =  0 

(1.3a) 

i>y(x.a)  =  GJ(x) 

,  0  <  x  <  » 

(1.3b) 

<!>(x  ,0)  =  Cx 

,  L  <  x  <  * 

(1.3c) 

i>y(x.O)  =  G‘(x) 

,  L  <;  x  c  ® 

(1.3d) 

where  the  constant  Cj  is  to  be  computed  by  iteration.  In  this  approximation, 
the  slip  lines  are  therefore  not  truly  free,  the  unknown  functions  G.!  being 
computed  from  the  solution  by  a  simple  differentiation.  The  remaining  boundary 

conditions  are 

<j>y(x»0)  =  f'(x)  = 

F(x)  ,  0  <  x  <  L 

(1.3e) 

$x(0,y)  =  h(>)  = 

H(Y)  . 

(1.3f) 
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Equation  (1.3f)  is  representative  of  the  initial  exit  velocity  profile  which 
conceivably  is  determined  by  the  upstream  duct  contour  and  stagnation  pressure. 

2.2  Far  Field 

To  complete  the  formulation  of  the  problem  for  subsonic  conditions  far 
downstream,  the  asymptotic  behavior  is  derived  in  this  section.  On  introducing  a 
Green's  function  G,  satisfying  homogeneous  Dirichlet  conditions  on  OA  and  CB  and 
homogeneous  Neumann  conditions  on  OC  and  AB  in  Fig.  2  with 

aG  =  6(P,Q)  =  6(x-£)s(y-n) 

where  Q(s,n)  is  the  source  point  and  P  is  the  field  point.  Green's  theorem 
applied  to  the  region  OQABCO  gives  the  following  integrodi fferential  equation 
for  <j> 

4 

*  =  l  Ii  . 

1  =  l 

where 

Ii  =  Sr  /  /  G(x»y;5.n)(3u2/3S)dn 

1  0  0 

Io  =  -I  ♦(5f0)  f|  (x,y;c,n)dc 

!3 =  _ci  [  (x’y;^0)d? 

a 

=  /  h(n)G(x,y;0,n)dn  • 

0 

The  quantity  G  may  be  obtained  from  the  cosine  transform  where 

oo 

5  =  /  G  cospx  dx 
0 

2  00 

G  =  -  /  5  cospx  dp  . 

11  0 
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The  subsidiary  equations  for  5  are 


d2S 

7 


n 2x  ./  \  x [  dG~l  /dG\  /dG\ 
p  G  *  6(y-n)cosp5  “>  -37  =  ( Hv  )  “  ( Hv ) 

l  •yJy=n  V  ,x/y=n+  \  vy=r 

■  cosp5  my-n  5  (GW  '  <s>y>„-  =  0 


5(p90;5,n)  *  6(p,a;e,n)  =  0  , 

implying  that 

G  a  ltnhPlaial.sl'lhpy„co_s£^ 

p  slnhpa  *  J  1 

-  sinhpn  sinhp(y-a)  cospg 

p  sinnpa  *  y  '  n  * 


(1.4a) 

(1,4b) 


Equations  (1.4)  can  be  inverted  by  a  treatment  of  appropriate  contour 
versions  for  the  inversion  integrals,  initially,  without  the  cosp£  factor  and, 
subsequently  including  it,  using  the  shift  theorem.  To  obtain  convergence  and 
exponential  decay  of  the  integrand,  the  appropriate  closure  for  the  contour  is  a 
large  semicircle  |p|  =  R,  R  +  <*>,  with  Im  p  \  0  for  x  \  0.  Summing  the  residues 
at  the  poles  p  =  niri,  n  =  (sgnx)(l,2,3,..,),  gives  the  following  final 
expression  for  G: 


-irG 


j,  n-1sin  noy  sin  non  {e~nct^x+^  +  e_Ra'x-sl } 


(1.5) 


a  =  Ti/a  . 


Equation  (1.5)  is  valid  for  y  ^  n»  ?»  and  can  be  summed  as  the  integral  of  a 
geometric  series  giving 


SC5055.21FR 


-2uG  =  S(a(x+£) ,a(n+y) )  +  S(a(x+£) ,a(n-y) ) 


-  S( otjx-5 1  »a(n+y))  -  S(a|x-?  |  ,a(n-y)) 


l-COSa(o-yj 

secha 

,x+c) 

X 

l-cosa(n-y)sechafx-e) 

l-cosa(n+y 

secha 

|x+d 

l-cosa(n+y)secha(x-5)J 

where 

S(A,B)  =  l  n-1e'nAcosnB  =  -£n | l-e~z | 
i 

=  -  j  *n(l-2e*AcosB  +  e'2A) 
z  =  A  +  iB  . 


(1.5') 


An  inspection  of  these  formulas  reveals  that  G  is  exponentially  small  as  x  +  », 
and  is  logarithmically  singular  at  the  source  point,  as  anticipated. 

Based  on  (1.5),  the  dominant  term  of  the  asymptotic  expansion  of  Ij  as 
x  +  oo  is  given  by  the  formula 

+  1  X  ^  ^  ? 

I,  4  ^  e"a  sincxy  /  sinandn  /  u  sinha?d^  (1.6) 

1  aK  0  0 

where  in  the  evaluation,*  the  contribution  of  the  simple  pole  of  G^  vanishes, 
and  integrals  of  the  form 

/do  /Va|x’^f(c,n)dC 
0  0 

and 

/3dn  /Va(x+5)fU,n)dS  (1.7) 

0  0 


*The  upper  limit  of  the  inner  integral  can  be  interpreted  as  ®  to  within  terms 
of  higher  order  involving  e_aX  as  x  +  ».  This  interpretation  is  made  in 
Eq.  (1.8). 
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arise.  The  multiplication  by  u2  of  the  asymptotic  expansion  represented  by 
(1.5)  as  x  +  m  and  Its  subsequent  Integration  with  respect  to  z  formally  gives  a 
development  dominated  by  these  Integrals.  Writing  the  Inner  Integral  of  the 
first  of  (1.7)  as 

^e-a|x-5lf(5#n)dc  -  e"aXjjXeaCf(5.n)d5  +  e°*/Ya?fU,n)ds  , 


and  if  the  e_ctX  factor  of  the  first  Integral  on  the  left-hand  side  is  indicative 
of  the  behavior  of  <j>  as  x  +  »,  then  u2  and  f  are  <^(e“2a?)  in  this  limit.  If  u2 
and  f  are  bounded  on  the  range  of  Integration,  the  first  Integral  converges  and 
the  second  Is  0(e_2otX)  as  x  +  <«>.  The  second  integral  in  (1.7)  requires  no  such 
decomposition  and  is  also  convergent;  hence,  (1.6)  follows.  Evaluating  the 
remaining  integrals  and  using  (1.5)  and  (1.5‘),  the  final  expression  for  the  far- 
field  is 


*  4  $pp  =  Cj (1-Y*)  +  Cppe'x*sinuY*  +  <?(e’2x*)  as  x  »  (1.8a) 


where 


cff  =  -Se  /  3ln*Y*dYJ  ^(?»Y)s1nhc*d5 
0  0 


®  o  1* 

-  —  /  H( Y) si n-jT Y*d Y  +  -Ar  j  »({•  ,0)cosh£*d?  -  2C.sinh  L*  (1.8b) 

*  0  0  1 


x*  =  hx/D/K  ,  Z *  =  ttC/D/K:  ,  L*  =  itL/D/K  ,  Y*  =  Y/D  . 

2.3  Numerical  Analysis 

The  numerical  procedure  Is  similar  to  one  first  developed  by 
Murman12  and  extended  by  Jameson^  and  Bailey  and  Ballhaus.*5  Briefly,  the 
transonic  potential  equation  in  divergence  form  is  discretized  by  using  central 
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differences  when  the  equation  is  elliptic,  and  backward  differences  when  it  is 
hyperbolic.  Thus,  we  may  write 

(K«x  *1)  +  (*Y)Y 

“  (l-ui)(K(*xH1/2"*xi-1/2)  “  *T  (<&xi+l/2"<txi-l/2) ^/pi 
+  wi-l(Kt*xi-l/2‘*x1-3/2)  '  ^  ($xi-l/2'$xi-3/2))/pi 


+  (*Yj+l/2  ‘  *Yj-l/2)/qj 


(2.1) 


-  (1-Ui) 


K 


/♦ur*i  +  VVij 

\ pi +1/2  pi -1/2  ) . 


,pi+l/2pi 


+  + 

Wl 

pi pi -1/2 . 


+ 


X 


/V*i-i 

\ Pi -1/2 


*  ¥  + 

♦r*i-i  _  *i-r*i-2 

pi-l/2pi  p i P i -3/2 


pi -3/2  / 

+  + 

,  *j+r*j 

qjqj*l/2 


+  + 

qjqj-l/2 


(2.2) 


where  the  missing  subscript  is  j  when  only  i's  are  present  and  i  when  only  j's 
are  present.  For  example,  <x> -  ^  j  and  $j+j  =  Also, 
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.  j? 


fir 


<  .* 


I  } 

< 


fT 

fi' 

K 


H 

If; 

fe 

fe-t. 


I? 


J 

p1-3/2  3  X1-l  “ 

* 

S- 

p1-l/2  =  X1  ‘  X 

i| 
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qj+l/2  =  Y1+l  ‘  Yj 
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VC 


3  K  -  Xti  fr*ll 1  +  $i  ~  , 
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0  if  the  point  (x^.Yj  is  elliptic 
1  1  J 

1  if  the  point  (x^)  is  hyperbolic  . 


Here,  the  iterations  are  viewed  as  steps  in  pseudotime  with  $  (NEW)  and  $  (OLD) 
values.  In  addition, 


!>  s  ( — J - 


pi-l/2pi  1  \p1p1-l/2  pipi-3/2/  ’\pip1-3/2, 


=3 

9£ 


g 
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These  definitions  guarantee  that  the  linearized  difference  algorithm  satisfies 
the  von  Neumann  stability  criterion.  (See  Jameson^  for  the  proof.) 

Overrelaxation  is  employed  in  the  elliptic  region  (y.  .  =  p.  .  .=  0) 

'  »  J  »  •  A  *  J 

only.  First,  define 


el  E  1/pipi-l/2  ’  e3  H  1/pi pi +1/2 


and 


e2  =  ei  +  e3 


Then,  the  elliptic  difference  expression  in  Eq.  (2.2)  given  by 


+ 

1+1  1 


+  + 

1  1-1 


pi+l/2pi  pi pi -1/2 
is  replaced  by 


=  e1*1-1  -  e2*i  +  e3$i+1 


el*i-l  "  *ie2/u)  '  $ie2(1‘1/(°)  +  e3*i+l 

where  u  is  the  overrelaxation  parameter,  i.e.,  1  <  w  <  2.  Note  that  if  oj  =  1, 
there  is  no  change  between  the  two  expressions. 

To  improve  stability  near  the  sonic  points,  especially  if  a 
discontinuous  wall  ooundary  condition  is  being  considered,  it  was  found 
necessary,  as  in  Bailey  and  Ballhaus,^  to  add  to  Eq.  (2.1)  the  term 


eAt 


Vi-i 


=  e 


K  -  ♦<)  -  (»i.i  -  «i.i ) 

<N  -  Vi>2 


where  e  is  chosen  to  be  in  the  range  0  <  e  <  .5. 
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Boundary  Conditions 


The  boundary  conditions  $(x,D)  =  0  and  $y(x,0)  =  f'(x)  may  be 
Incorporated  into  the  numerical  scheme  by  using  the  same  techniques  described 
In  Murman  and  Cole.1^  Discussed  first  will  be  the  boundary  condition  at  the 
jet  exit  x  =  0,  which  may  be  one  of  two  types:  (A)  subsonic  at  the  jet  exit, 
or  (B)  partially  or  completely  supersonic  at  the  jet  exit.  The  far  field 
boundary  condition  will  be  treated  later. 


(A)  Subsonic  at  the  Jet  Exit 

Here  the  boundary  condition  is 


*X(0,Y)  =  H(Y) 


0  <  Y  <  D 


where 


K  -  (y+l)H(Y)  >0  for  0  <  Y  <  D  . 

In  this  case,  let  x^2  =  0,  -  x^2  =  Ax/2,  Mg  =  0,  and  y-j  =  0; 
it  is  also  required  that  x2  -  x-j  =  Ax.  Then,  the  derivatives  in  Eqs.  (2.1) 
and  (2.2)  become 


[K^x3/2  ‘  ^x  1/2 ^  ‘  V"  ($x3/2  '  *xl/2 )]/AX 


[k  (®x3/2  +  *xl/2^][®x3/2  '  ®xl/2 ] / 

*  H<¥))][~Tr '  H(v)]/S: 
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(B)  Partially  or  Completely  Supersonic  at  the  Jet  Exit 

For  this  case,  two  boundary  conditions  are  required  at  x  =  0;  namely, 

*X(0,Y)  =  H( Y)  0  <  Y  <  D 


and 


4(0, Y)  =  g(Y)  for  Y  e  [0,0] 


where 


K  -  (y+l)H(Y)  <  0  . 

For  all  points  Y  e  [0,0]  in  which  K  -  (y+l)H(Y)  >  0,  g(Y)  need  not 

exi st. 

If  a  point  (x^,Yj)  is  elliptic,  the  discretization  given  in  case  (A) 

is  used.  On  the  other  hand,  if  (x-,,Y.)  is  hyperbolic,  it  is  assumed  that  the 

J 

grid  may  be  extended  to  the  left  by  Ax/2,  and  that  Xq  =  0,  x =  -Ax/2,  and 

V,  ■  *<°.V  -  g<vj). 

Using  Taylor's  tneorem, 

‘-Id  =‘0j  -  M*x(0.»j)/2  ■  9(Vj)  -  AxH(rj)/2  . 

These  values  for  and  4_i>j  <nay  now  be  substituted  into  Eq.  (2.2)  in  the 
normal  way,  and  line  relaxation  may  be  applied  to  the  first  column  of  unknowns 
along  x  =  xj. 

The  far  field  boundary  conditions  given  by  Eq.  (1.8)  contain  two 
unknown  constants,  Cj  and  Cpp,  which  must  be  determined  in  an  iterative 
fashion.  The  basic  technique  holds  Cpp  fixed  while  changes  until  the 
solution  converges.  Then,  Cpp  is  updated  by  evaluating  the  integrals  in 
Eq.  (1.8b),  and  the  procedure  is  repeated  until  Cpp  also  converges.  In  order  to 
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determine  C^,  the  mesh  network  is  swept  from  left  to  right  by  using  line 
relaxation.  After  the  potential  is  computed  on  the  line  x  =  L,  extrapolation 
of  the  interior  points  yields  $(1,0)  s  $  If  it  is  assumed  only  that  $  is 
continuous  at  x  =  L,  then  =  $(L,0)  =  4>(x,0)  for  x  >  L,  which  guarantees 
that  $x(x,0)  =  0  for  x  >  L. 

Singular  Behavior  Near  Mall  Trailing  Edge 

It  is  physically  plausible  that  a  Kutta  condition  given  by 

♦x(L-.°)  =  <t>x(L+»0)  -  0  (3) 

is  satisfied  by  the  solution  for  trailing  edge  neighborhoods  In  unmixed  flow. 
Because  of  (3),  the  nonlinear  term  in  (1.2)  can  be  assumed  negligible,  and  4,  is 
locally  harmonic  in  the  scaled  variables.  Let 

z  =  x*  +  iy  ,  x*  =  x  -  L 

e  =  arg  z  ,  r  =  mod  z 

w(z)  =  u(x*,y)  -  1v(x*,y)  =  complex  velocity 

where 

u  “  V  .  V  =  4>y  • 

If  the  boundary  conditions  are  locally  linearized  near  the  point  z  =  0, 
then 

v(x*,0)  F(L)  so)  ,  x*  i  0  (4a) 

u(x*,0)  =0  ,  x*  >  0  .  (4b) 
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To  dominant  order,  a  sufficient  condition  to  satisfy  (3)  and  (4)  near  the  origin 
is  that  w  has  the  following  branch  point  behavior 

w  =  i (to  +  B/z)  as  z  >  0  ,  (0  <  9  <  u)  (5) 

where  B  is  a  real  constant  to  be  determined  by  matching  with  the  outer  nonlinear 
solutions.  Equation  (5)  implies  that 

4>  -  <(>L  i  -wy  -  Br3^2sin3e/2  .  (6) 

Several  examples  to  be  discussed  indicate  that  the  approach  described  previously, 
in  which  $  is  maintained  continuous  at  the  wall  trailing  edge,  gives  numerical 
solutions  that  satisfy  the  Kutta  condition  (3).  However,  a  rigorous  proof  that 
this  is  an  implication  of  the  algorithm  has  not  been  attempted.  A  similar 
procedure  has  been  used  by  Krupp17  to  satisfy  the  Kutta  condition  in  the 
solution  of  the  transonic  small  disturbance  lifting  airfoil  problem. 

2.4  Results  and  Discussion 

In  addition  to  the  assumptions  given  in  the  Introduction,  the  analysis 
previously  described  is  not  directly  applicable  to  choked  flows  where  upstream 
and  downstream  conditions  are  decoupled.  Sonic  zones  comprising  the  entire 
vertical  dimension  of  the  flow  field  are  thereby  excluded.  However,  the 
foregoing  methods  can  be  extended  to  handle  such  cases. 

A  number  of  examples  will  now  be  considered.  For  these  cases,  the 
associated  wall  displacement  functions  and  exit  velocity  distributions  H(Y)  are 
given  in  Table  1.  For  most  of  these  cases,  K  is  unity  and  D  will  take  on  this 
value  for  the  remainder  of  this  paper. 

Wall  pressure  distributions  for  the  convex  ramp  comprising  Case  1  are  shown 
in  Fig.  3  for  K  =  1  and  K  -  1.46.  By  interpreting  these  results  as  distributions  for 
different  final  M  ‘s  downstream  but  with  the  same  5,  the  decrease  in  to  leads  to 
upstream  motion  of  the  terminating  shock  but  a  preservation  of  the  shock  strength. 
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Table  1.  Wall  Jet  Cases 


Case 

f  =  Wall  Shape  Function 

H(Y) 

Remarks 

1 

f  *  fj  =  0  .  0  <  x  <  1 

=  -(1-x)2  ,  1  <  x  <  2 

*  3-2x  ,  2  <  x  <  L 

.075 

0  <  Y  <  1 

L  =  3.8  for  this  and  all 
other  cases 

2 

f  »  f2  *  -fl 

3 

f  *  f3  =  0  ,  0  <  x  <  1 

=  [x2-l+2(l+l) ( l-x)]/2L  ,  1  <  x  <  L 

Has  discontinuous  slope 
at  x  =  1 

4 

f=fjj=0  ,  0  <  x  <  1 

.  Mtdi  cos  llftl  -  1  ,  1  <  x  <  L 

Has  reflex  curvature  on 
curved  ramp  portion 

5 

f  *  f5  =  -x2  ,  0  <  x  <  i 

«  l-2x  ,  1  <  x  <  2 

=  -3  ,  2  <  x  <  L 

H  =  1  0  <  Y  <  1/2 

■=  0  ,  1/2  <  Y  <  1 

f  =  0.2  on  0  <  V  <  1/2 

25 

C3471A/J  bs 


r«<w> WOUU j Si'W. W. • . »%  A  !»«&«>  •■-  — 


* 


Rockwell  International 

Sc  tone*  C#nt*r 


SC5055.21FR 

There  is  a  smooth  acceleration  to  critical  conditions  with  the  location  of  the 
sonic  line  established  at  the  curvature  discontinuity,  x  =  1.  These  calcula¬ 
tions,  which  are  typical  of  the  other  cases,  cost  approximately  $60  on  the 
Berkeley  7600  and  ran  15-30  CP  seconds.  Figure  4  shows  a  close-up  of  the 
pressures  near  the  trailing  edge.  The  dashed  line  has  a  slope  proportional  to 
/**,  appropriate  to  the  singular  behavior  given  by  Eqs.  (5)  and  (6)  and  the  Kutta 
condition  (3).  Isobars  shown  in  Fig.  5  are  consistent  with  these  remarks  and 
demonstrate  the  satisfaction  of  homogeneous  pressure  boundary  conditions  on  the 
slip  lines.  Because  of  the  weakness  of  the  singularity,  e.g.,  <j>xx  -  r"*^  as 
r  -v  0,  special  numerical  treatments  such  as  those  of  Woods*8  were  not  used. 

In  Figs.  6  and  7,  rapid  decay  of  tne  disturbances  Is  indicated.  The 
relaxation  length  for  this  decay  from  (1.8b)  with  D  =  K  =  1,  is  n.  This 
exponential  decay  is  typical  of  flows  confined  by  jet  boundaries  and  is  much 
more  potent  than  for  bodies  In  unbounded  fields.  Qualitatively  similar  effects 
have  been  discussed  by  Murman,19  and  Pinzola  and  to20  in  connection  with  tunnel 
wall  interference  on  transonic  airfoils.  The  distinction  between  confined  and 
unconfined  flows  can  be  appreciated  by  an  interpretation  of  the  exponential 
series  due  to  (1.5)  arising  in  the  far  field  developments  dominated  by  (1.8). 
Because  of  linearity  of  the  far  field,  this  series  is  directly  related  to 
expansions  occurring  in  analogous  incompressible  flow  problems,  in  which  a 
singularity  is  reflected  between  free  pressure  boundaries.  This  yields  an  image 
development  in  which  the  strengths  alternate  in  sign  to  satisfy  the  slip  line 
condition.*  Thus,  the  relaxation  to  uniform  conditions  downstream,  which  must  be 
consistent  with  homogeneous  conditions  on  the  slip  lines,  produces  a  more  rapid 
decay  than  found  in  unconfined  flows. 


*Such  a  series  can  be  summed  by  recognizing  that  it  is  a  partial  fraction  ex¬ 
pansion  of  a  hyperbolic  function  which  is  exponentially  small  as  x  +  ®,  con¬ 
sistent  with  developments  such  as  (1.8a). 
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In  Figs.  8a  and  8b,  the  upper  and  lower  slip  lines  obtained  from  inte¬ 
gration  of  (1.3b)  and  (1.3d)  are  given  for  Case  1.  It  is  evident  from  Fig.  8a 

that  the  curved  surface  in  this  approximation  turns  the  flow  so  that  the  streams 
are  parallel  for  x-*®.  It  is  evident  from  Fig.  8b  that  in  the  near  field,  this 
is  not  quite  the  case.  Asymptotic  parallelism  can  be  established  for  subsonic 
conditions  far  downstream  by  integration  of  the  small  disturbance  continuity 
equation.  Thus 

S(x)  "  G2(x)  =  4  MKu(x*Y)  '  U^dY  (7) 

and,  since  u=  £?(e‘x*),  this  expression  shows  that  Gj(x)  -*•  G^(x)  as  x 

Equation  (7)  can  also  be  obtained  by  differentiating  (1.8a)  with  respect  to  Y 
and  using  (1.3b)  and  (1.3d).  For  the  case  of  a  free  jet  with  a  symmetrical  exit 
Mach  number  profile  function  H,  i.e.,  H( Y-l/2)  =  H(l/2-Y),  Gj  =  -G2  and  the 
divergence  theorem  or  integration  of  (7)  between  x  =  0  and  »  gives  the  displace¬ 
ment  of  the  jet  infinitely  downstream  as 


where  G^(0)  =  0  has  been  used. 

Insight  into  the  mechanisms  causing  acceleration  to  supercriticality 
can  be  obtained  from  the  flow  direction  field  for  Case  1  shown  in  Fig.  9.  For 
clarity,  all  isocline  slopes  have  been  magnified  by  a  factor  of  100,  and  only 
the  entrance  section  0  <  x  <  1  is  depicted.  The  expansion  around  the  curved 
ramp  on  the  interval  1  <  x  <  2  leads  to  upstream  influence  in  the  subsonic 
region  which  turns  most  of  the  flow  downward  upstream  of  x  =  1,  producing 
throats  and  acceleration  due  to  the  stream  tube  contraction  required  by  the  zero 
slope  boundary  condition  in  that  region.  Also  indicated  is  the  "ballooning"  due 
to  the  singularity  occurring  at  (0,1),  the  top  point  of  the  jet  exit  station. 

In  contradistinction  to  the  trailing  edge  where  u  0,  local  linearization 
cannot  be  used  to  characterize  the  flow  behavior  in  this  region,  and  some  local 
similarity  solution  must  be  sought,  presumably  of  the  foi  n  <*,  =  xuf(y/x^)  where 
a  and  3  are  exponents  to  be  determined. 
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By  contrast,  the  concave  shape  shown  in  Fig.  10  produces  the  antici¬ 
pated  compressive  deceleration,  which  is  also  depicted  in  Fig.  11. 

The  effect  of  a  slope  discontinuity  is  indicated  in  Fig.  12.  It  is 
evident  that  the  numerical  method  accurately  locates  the  initiation  of  the  sonic 
line  at  the  point  (1,0)  where  the  flow  is  "tripped"  to  criticality  by  the 
acceleration  singularity  at  this  location.  In  most  other  respects,  the  pressure 
distribution  is  similar  to  that  for  Case  1. 

In  Fig.  13,  the  effect  of  a  reflex  curvature  in  accelerating  the  recom¬ 
pression  process  is  shown.  As  related  to  comparable  turning  and  wall  deflection 
treated  in  Case  1,  the  strength  of  the  terminating  shock  is  considerably  in¬ 
creased,  as  is  the  magnitude  of  the  pressures  near  the  trailing  edge. 

The  effect  of  mixed  flow  conditions  at  the  exit  is  shown  in  Figs.  14-16. 
Here,  the  function  H  as  well  as  <t(0,Y)  comprise  the  Cauchy  data  needed  to 
properly  pose  the  hyperbolic  portion  of  the  initial  manifold.  Since  the  verti¬ 
cal  velocity  4>y(0,Y)  can  be  derived  as  a  tangential  differentiation,  the  Cauchy 
data  connotes  specification  of  the  additional  velocity  component  for  supersonic 
portions  of  the  jet  exit  station.  Figure  16  indicates  that  in  addition  to  the 
usual  terminating  shock,  the  transition  from  hyperbolic  to  elliptic  flow  occurs 
across  a  weak  shock  emanating  from  the  specified  u  discontinuity  at  (0,1/2). 

In  Fig.  17,  the  behavior  of  centerline  pressures  for  various  free  jet 
cases  is  shown.  The  monotone  subcritical  behavior  exhibited  by  these  nonlinear 
cases  has  not  yet  been  corroborated  by  rigorous  proof  based  on  the  boundary 
value  problem  with  subsonic  exit  and  downstream  conditions.  For  lir,eariz_d 
subsonic  flow,  this  property  is  obvious  from  the  maximum  modulus  theorem.  It 
should  be  noted  that  in  the  free  jet  problem  discussed  here,  specified  mass 
flow,  pressure  ratio,  and  final  Mach  number  M^  uniquely  determine  6,  the  scale 
parameter  for  the  jet  displacement. 

For  a  validation,  a  comparison  of  numerical  and  Prandtl-Glauert  free  jet 
solutions  for  u(x,l/2)  for  H  =  0.35  and  K  =  10  is  shown  in  Fig.  18,  where  the 
analytical  solution  obtained  either  from  summed  eigenfunction  expansions  or 
transforms  is: 
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Fig.  10  Pressures  along  wall  for  Case  2,  fo  =  normalized  wall  ordinate 


Fig.  12  Pressures  along  wall  for  Case  3,  f2  =  normalized  wall  ordinate 
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(9a) 


(9b) 


and 


u  .  IS  tan-l  sin/ 
sinhx1 


_  H/K  „„  coshx*  +  cosy* 

*  ”  Xr'l  1*“  .1. 

coshxT  -  cosy1 


x*  =  tix//K  ,  y*  =  ttY 


The  slight  discrepancies  shown  in  Fig,  18  presumably  derive  from  the  small 
nonlinear  effect  associated  with  the  finite  K  value  and  from  truncation  errors  of 
the  discretizations  which  are  only  approximately  second  order  for  a  non-uniform 
elliptic  mesh.  The  associated  universal  slip  line  curve  is  obtained  from  the 
following  integral  of  (9b) 


G4-  .  >-  -  *  2  n.tACtntllltl!/1.!!  .  X  ctnh-hcoshu)  (10) 

HK  2  it  0  (2n  +  l  f  2 


where  the  daggers  have  been  dropped.  Equation  (10)  has  the  following  limiting 
behavior 


^1  (x)  p  (  /I  \  2  ) 

-^lj7 —  =  ~2  | -x jinx  +  yjr  2o4+lj  x  +  f'(x  )> 


^i(x)  i  4  _x  _2x 

“HIT-  =  2  "  ~2  6  +  0  (xe  ^ 


as  x  +  0 


as  x  +  » 


and  is  plotted  in  Fig.  19.  The  asymptotic  half  width  thus  checks  that  given  by 
(8)  when  the  nonlinear  second  tenn  in  the  integral  is  omitted. 
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3.0  I NV I  SC  ID  MODEL  FOR  UPPER  SURFACE  BLOWING 

In  the  previous  section,  a  submerged  jet  problem  related  to  tangential 
blowing  was  studied.  Actually,  the  jet  structure  investigated  in  this  connec¬ 
tion  can  be  realized  only  approximately  for  large  values  of  the  jet  momentum 
flux.  For  more  moderate  values,  it  is  necessary  to  account  for  interaction  with 
the  external  stream.  This  implies  that  a  coflowing  rather  than  submerged  jet 
must  be  considered  adjacent  to  the  wall  and  that  the  associated  perturbations  of 
flow  external  to  the  jet  be  treated. 

In  Fig.  20a,  a  schematic  of  an  upper  surface  blown  airfoil  is  shov/n 
embedded  in  a  transonic  flow,  with  a  supersonic  wall  jet  tangentially  injected 
along  its  upper  surface.  The  inviscid  problem  for  this  configuration  is 
formulated  in  Fig.  20b,  based  on  a  heuristic  generalization  of  the  jet  flap  model 
of  Ref.  21.  Conditions  across  the  jet  and  jet-wake  are  not  indicated.  Relevant 
to  this  formulation,  we  will  discuss  certain  aspects  of  the  coflowing  jet  flow 
which  are  pertinent  to  establishing  the  appropriate  boundary  conditions  for  the 
external  flow.  We  will  also  treat  other  properties  that  have  significance  for 
non-aeronautical  applications  in  addition  to  the  present  one. 

Considering  the  coflowing  wall  jet,  we  divide  its  consideration  into 
regions  such  as  R?,  Rj£,  and  domains  outside  of  these  neighborhoods  on  and  off 
the  airfoil.  The  latter  domains  are  inviscid  boundary  layers  with  a  structure 
consistent  with  a  balance  solely  between  pressure  and  centrifugal  forces.  A 
systematic  asymptotic  theory  for  the  fine  structure  of  these  layers  is  given  in 
Section  3.2.1.  In  other  sections  and  Refs.  5-7,  properties  of  Re  and  R^  are 
discussed.  For  these  regions,  the  effects  of  streamwise  acceleration  and  wave 
interaction  phenomena  controlled  b/  a  length  scale  of  the  order  of  the  jet 
thickness  significantly  alter  this  balance.  To  shed  light  on  the  transition  of 
these  axial  boundary  layers  to  the  other  zones,  Section  3.1  treats  an  “inner" 
problem  associated  with  the  coflowing  wall  jet  relevant  to  the  matching  of  the 
regions  Re  and  to  the  jet,  wake,  and  external  flow.  Section  3.1.1  also 
discusses  selection  rules  to  define  the  nature  of  the  mixed  flow  structures  that 
can  arise  in  the  wall  jet.  Numerical  studies  illustrating  these  structures 
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which  produce  significant  spikes  observed  In  blown  transonic  airfoil  pressure 
distributions  are  provided  subsequently  In  Section  3.1.3.  In  Section  3.2.2, 
other  aspects  are  discussed  relevant  to  and  the  formulation  for  the  global 
description  for  the  blown  airfoil  problem.  Results  are  given  for  typical  cases 
in  later  portions  of  this  report. 


3.1  Coflowirig  Transonic  Mall  Jets  -  An  "Inner11  Problem 

In  Fig.  21a,  the  appropriate  formulation  of  a  problem  In  which  the  re¬ 
gions  Re  and  RyE  are  merged  In  the  tangentially  blown  jet  region  of  Fig.  20a  Is 
shown.  As  Indicated  previously,  we  are  concerned  with  the  transition  of  these 
regions  in  which  streamwise  gradients  are  Important  to  other  regions  which  are 
dominated  by  gradients  normal  to  the  jet.  For  the  former,  we  note  that  for 
supersonic  and  transonic  portions  near  the  jet  exit,  the  characteristic  longitu¬ 
dinal  scale  is  of  the  order  of  the  wavelength  of  the  reflected  wave  pattern  from 
the  slip  lines.  To  treat  the  longitudinal  scales  properly,  a  KG  limit,  which  in 
the  notation  of  Fig.  20  has  (x^E  -  xe)/64^3,  x/S4^3,  [y  -  6b^(x)]/6, 

K  =  (1  -  R“:)/6^3,  K  =  (1  -  ) / 6 2 ^ 3  fixed  as  6+0,  is  utilized.  Here,  is  the 

G  G  oo  G 

average  of  the  square  of  the  exit  Mach  number  distribution,  and  is  the  freestream 
value  of  the  Mach  number.  In  this  formulation,  a  new  excess  pressure  ratio 
parameter  arises  naturally,  which  is  given  by 


U)  = 


,  fixed  in  the  KG  limit 


where  is  the  freestream  ambient  static  pressure  and  Pe  is  the  corresponding 
static  value  at  the  jet  exit.  Essentially,  the  approximation  of  slightly 
different  velocities  between  the  jet  and  external  flow  is  used  to  study 
important  flow  features.  In  practical  applications  at  transonic  speeds,  this 
could  be  realized  with  nearly  choked  blowing  nozzles. 

The  formulation  shown  in  Fig.  21a  can  be  derived  in  the  usual  way  by 
substituting  the  appropriate  asymptotic  expansions  for  the  velocity  potential 
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Fig.  21b  Jet  region  -  thin  membrane  theory. 
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Into  the  full  potential,  and  energy  equations  and  the  boundary  conditions.  On 
the  boundaries  of  the  cowl  B'V*  and  CV,  the  latter  Involve  tangency  or  x.  .& 
flow.  On  the  slip  lines  C'U'  and  CU,  they  Involve  continuity  of  pressure  and 
vertical  velocity.  At  the  jet  exit,  Oirichlet  or  Neumann  conditions  are 
specified  if  the  exit  is  subsonic.  If  the  exit  is  supersonic,  Cauchy  data 
are  required.  This  model  provides  a  means  for  studying  the  coupling  of  the 
external  and  internal  flows  and  for  understanding  the  role  of  this  interaction 
on  the  development  of  the  flow  on  the  rear  portion  of  the  airfoil.  Extension  to 
axially  symmetric  cases  also  appears  feasible.  Intrinsic  in  the  formulation  is 
the  representation  of  the  internal  and  external  flows  by  separate  perturbation 
potentials  $  and  «,  respectively.  The  appropriate  asymptotic  expansions  for  each 
are : 

JET 

~  $  •;  2 

if-  =  x  +  6  4>(x,y)  ►  ••• 
e 

4>i  =  exact  jet  potential 

x  =  x/64/3,  y  =  (y  -  «ba(x))/«,  Ke 
fixed  as  6+0 

Ue  =  average  velocity  at  jet  exit 
BLENDING  LAYERS  (IMMEDIATELY  ABOVE  AND  BELOW  JET) 

♦  d  2 

jp  =  X  +  6  4>(x,y)  +  •••  (3) 

x,  y,  K  =  (1  -  M^)/6^3,  a)  fixed  as  5+0 
00 

U  "  freestream  velocity  in  external  flow 


(1) 

=  (i  -  Rg)/62/3,  o)  =  (r  ■  (2) 

00 
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EXTERNAL  FLOW 

-{p-  =  x  +  62/34>(x,y)  +  ...  (4) 

1/3 

x,y  =  6  '  y,  K,  to  fixed  as  6-.0  . 

If  *  represents  the  velocity  potential  in  any  of  the  three  regions  delineated 
above,  the  velocity  is  given  by  74.  On  the  cowl  surfaces  and  the  slip  lines 
denoted  generically  as  S(x,y)  =  0,  we  have 


q  •  VS  =  0 


as  the  tangency  boundary  condition.  In  connection  with  the  slip  lines,  note 

that  6b  (x)  =  (7).  Moreover,  the  perturbation  pressures  for  the  jet  and 

i  J 

blending  layer  are,  respectively: 


JET 

pi  =  1  -  Y62/3P i (x ,y)  +  ... 
e  J 

where  Pe  is  the  mean  pressure  at  the  jet  exit. 


BLENDING  LAYER 

•p  =  1  -  y62/3pB(x,y)  +  ••• 

00 

where  is  the  ambient  pressure  in  the  freestream. 


EXTERNAL  FLOW 


l  -  Y$2/3p(x,y) 


+  •  •  • 


(5a) 


(5b) 


(5c) 


52 

C347 1 A/ j  bs 


Rockwell  International 

Science  Center 


SC5055.21FR 

By  virtue  of  pressure  compatibility  on  the  slip  lines,  the  relation 


(6) 


holds  along  the  slip  lines  8' U '  and  CU  in  Fig.  21a,  which  upon  tangential  inte¬ 
gration  yields  the  relations  indicated.  The  constants  Cu  and  are  similar  to 
circulations  and  are  obtained  from  the  difference  of  the  potentials  at  the 
trailing  edges  B'  and  C.  The  trailing  edge  behavior  and  jet  turning  angle  are 
obtained  from  (6)  and  the  compatibility  relations  <P-  =  <)>-  on  the  slip  lines. 
Local  nonuniformities  occur  at  the  trailing  edges  which  are  anticipated  as  weak 
singularities.  "Inner  inner"  expansions  for  these  would  require  the  solution  of 
the  full  potential  equations  in  these  neighborhoods. 


To  complete  the  specification  of  the  problem,  a  far  field  <j>pp  is  neces 
sary.  This  function  should  have  the  property  that  it  asymptotically  matches 
with  the  outer  flow  on  the  boundary  of  the  region  R0.  A  procedure  which  pro¬ 
vides  the  necessary  nonlinear  feedbacks  with  the  near  field  of  the  problem  in 
Fig.  21a  involves  the  asymptotic  solution  of  a  relevant  integrodi fferential 
equation.  A  heuristic  far  field  derivation  in  anticipation  of  this  more  satis¬ 
fying  approach  is  described  in  the  next  section. 

MATCHING  OF  EXIT/TRAILING  EDGE  REGION  WITH  EXTERNAL  FLOW 


Because  of  the  anticipated  jump  in  the  vertical  velocity  at  the  trail¬ 
ing  edge  associated  with  a  jet  flap  or  USB,  we  assert  that  the  inner  limit  of 
the  external  flow  solution  which  is  assumed  to  be  locally  harmonic  in  strained 
variables  in  that  neighborhood  has  a  logarithmic  singularity  in  the  complex  per¬ 
turbation  velocity.  To  simplify  the  discussion,  we  consider  matching  in  the 
*  co- text  of  a  jet  flap  (xe  =  xTp)  for  the  present.  Generalization  to  more  com¬ 

plicated  USB  cases  should  present  no  major  new  difficulties.  As  a  consequence 
f  of  the  logarithmic  singularity  assertion,  the  complex  potential  behaves  like 

z  log  z,  as  z+0,  where  z  =  x  +  i/K  y  in  the  vicinity  of  the  trailing  edges  as 

i 

M 
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they  appear  in  the  external  formulation.*  Accordingly,  in  the  notation  of  (6), 
we  assume  that 


<j>  =  A{x  log  r  -  A  ya}  +  •  (7) 

o  p  o  i  jiT  rt 

where  A  is  a  constant,  r  -  x  +  Ky  ,  and  9  =  tan"  ~ -y-  .  In  a  corresponding 
manner,  we  assert  that  the  outer  limit  of  the  blending  layer  solution  for  $ 
associated  with  the  expansion  (3)  is  given  by 

<j>  =  a{x  log  r  -  /i<  y  0}  +  ...  (3) 


where  a  is  another  constant. 

For  matching,  we  introduce  intermediate  variables  x^  and  y^  given  by: 

VlRT  *  yn  =  ^7  fixed  as  (9) 

with 

6  4/3  «  K  «  1 
6  «  n  «  6'1/3 

where  the  «  sign  signifies  "order  of." 

In  the  intermediate  limit  (9)  which  implies  that 

x  =  r,xr  ->  0 

y  =  o1/3nyn  +  0 


*For  near  critical  trailing  edges  this  assumption  is  invalid. 
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x’  =  7h 


+  CO 


y  n 

y  =  — —  +  “ 
J  5 


the  variables  and  various  terms  appearing  in  (7)  and  (8)  have  the  following 


appearance 


-1  /y  -  -  /K  ny  6 

9  =  tan"1  1  +  tan-1 - -- 

x  X 


i  /F  w  i  /K  ny  6 

0  =  tan-1  >  tan"1  — — n-- 
x 


/  52x2  n2y2 

log  F  - log  v  -m +  -/ 


log  r  =  log  J X x^  +  K52/3n2y2  • 


(10a) 


(10b) 


(10c) 


(lOd) 


On  the  basis  of  a  comparison  of  the  various  elements  in  Eqs,  (10),  we  find  by 
virtue  of  (7)  and  (8)  that  the  perturbation  potentials  in  b^th  regions  compare 
roughly  as  follows,  where  the  <=>  indicates  asymptotic  correspondence: 

6^3<j>  <=>  52<Ji  +  4$  0 o?  *)  x  (11) 


providing  the  limit  is  taken  as  — fixed  =  constant  =  x  as  5+0.  More 
precisely,  if  the  external  and  inner  flow  solutions  are  written  in  terms  of 
intermediate  variables,  we  obtain  the  following  correspondence: 
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As*7'5}™  log  +Tx7c7y2’  -  /K  xcy  9} 

s  s  n  n 


rx  1 

<=>  5?a  ~Kft  log  ^2xp“x^Ty?  -  log  5  -  /K  -J^y^ 

5  '6 


+  4(«2">9«)^j  . 
6 


The  last  term  is  a  compensation  term  added  to  achieve  matching.  Presumably,  an 
intermediate  solution  could  be  developed  to  eliminate  this  rather  foreign  look¬ 
ing  entity  whose  physical  significance  is  not  well  understood  at  this  time. 

Aside  from  this  rather  minor  difficulty,  the  result  contained  in  (12)  demon¬ 
strates  that  a  numerical  solution  of  the  problem  in  Fig.  21a  will  match  with  the 
solution  of  a  global  jet  flap  problem  near  the  trailing  edge  providing  A  =  a  and 
the  far  field  for  the  former  is  given  in  accord  with  (12). 


"THIN  DOMAIN"  SPENCF.  REGION 

Another  matching  is  required  downstream  of  the  trailing  edge,  involving 
the  confluence  of  the  wavelike  region  and  the  "Spence"  zone  where  pressure  and 
centrifugal  forces  equilibrate.  This  regime  is  reached  in  the  limit  of  x,y 
fixed  as  and  leads  to  the  "thin  domain"  problen  shown  in  Fig.  21b.  Here, 
the  streamwise  derivatives  are  negligible  in  the  dominant  approximation. 

Matching  with  the  external  flow  takes  place  by  "patching"  across  the  slip  line 
with  the  boundary  conditions  of  compatible  normal  velocity  and  pressure  as  shown 
in  the  figure.  The  inner  limit  of  the  thin  layer  solution  represents  the  jet 
far  field  of  the  problem  associated  with  Fig.  21a. 

3.1.1  A  Basic  Transonic  Coflow  Problem 

Many  of  the  fundamental  ideas  associated  with  the  problem  depicted  in 
Fig.  21b  are  contained  in  the  unit  problem  shown  in  Fig.  22.  Here,  a  jet  of  ve¬ 
locity  Ue  exhausts  into  a  coflowing  stream  with  a  freestream  velocity  U.  A 
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major  simplification  over  the  problem  given  in  Fig.  21a  is  now  that  the  cowl  PBO 
is  a  flat  plate  internally  and  externally  as  well  as  the  bottom  wall  surface  LSR. 
Moreover,  the  lower  cowl  boundary  VC  in  Fig.  21a  is  replaced  by  an  infinite  wall 
y  =  -1  in  Fig.  22  for  this  idealization.  We  assume  that  conditions  are  known  at 
some  exit  station  AB,  which  could  be  at  x  *  -«  in  the  coordinate  system  indicated, 
and  could  correspond  to  a  settling  chamber  location.  Note  in  this  connection 
that  a  far  field  analysis  indicates  that  as  station  AB  tends  to  x  =  the 
assumed  Neumann  condition  is  asymptotically  equivalent  to  a  homogeneous 
Dirichlet  one.  In  particular,  with  the  problem  for  $  based  on  Fig.  22,  a  far 
field  representation,  obtained  from  an  iterative  solution  procedure  regarding  the 
nonlinear  term  of  the  Karman  Guderley  equation  as  a  weak  forcing  term,  is 


irx//R7  A.2  2ttX//K~ 

9(x, y)  =  Aj  e  cos  *y  +  g±-  (y+l)x  e  cos  2*y 


37? 


8  K 


2«x//K 

e  e  +  A2 


2ttx//K 

e  cos  2ny  + 


•  •  • 


as  x  -*■  -<*>  . 


Defining  the  average  of  a  flow  quantity  g  at  the  station  x  =  x0  as 

0 

<9(x0,y)>  =  avera9e  of  g  at  x  =  xQ  =  /^>(x0,y)  dy 

we  have  in  particular,  using  the  foregoing  asymptotic  solution  specialized  to 
x0  =  1. 

9  /  ii\  2  -2tt//kT 

<y(-i.y)>  •  a.2  e  . 

1  8 

Moreover,  we  find  that 


<?x(-l.y)>  s  7!1  <?{-i.y)> 
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An  interpretation  of  the  constant  A  appearing  in  the  foregoing 
expressions  is  obtained  by  using  Green's  theorem  to  derive  a  nonlinear  integro- 
differential  equation  for  <p .  From  the  formulation  of  Fig.  22,  it  is  possible  to 
evaluate  the  boundary  terms  on  the  boundaries  of  the  rectangular  domain  AB  and 
PR.  The  appropriate  integral  theoren  is 


0  «  #  j°  U2M  -  dn 

^  -1  ^  o  -®  F-  ^ 


-  -  /  b'(c)  in  sinh  d^  +  i  x  , 


where 


Approximation  of  this  expression  for  x  +  ±  <*>,  with  interrelationships  betv.-een 
b(a>)  and  c  to  be  derived  subsequently,  is  consistent  with  the  asymptotic 
expansion  of  <p  for  x  -*•  -*>  previously  given  and  the  result 

c 

<j>  =  j  x  as  x  i  oo 

where  the  partial  integration  boundary  term  in  the  first  integral  of  the  Green's 
theorem  expression  dominates  the  double  integral  and  is  approximately  +  cxc/2  as 
x  •>  »,  and  the  contribution  of  the  lower  limit  vanishes.  This  change  in  signs 
is  of  decisive  importance  in  obtaining  the  desired  upstream  and  downstream 
boundary  conditions.  Furthermore,  the  Green's  theorem  expression  shows  that 


A1  =  /  b'U)  e"715  dc  , 


i.e.,  A^  is  a  moment  of  the  source  distribution  associated  witn  the  deflected 
slip  line. 
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Returning  to  the  configuration  at  hand,  we  note  that  in  some  sense, 
this  system  is  a  simple  idealization  of  an  inviscid  ejector.  A  sketch  of  a 
wind  tunnel  realization  is  also  shown  in  Fig.  22.  An  analytical /computational 
solution  has  been  developed  to  study  this  problem,  with  the  objective  of 
addressing  the  following  issues  associated  with  the  problem  of  Fig.  21b: 

1.  Under  what  circumstances  is  the  jet  far  field  supersonic  if  the 
initial  conditions  are  subsonic? 

2.  What  is  the  degree  of  downstream  penetration  of  an  initially 

supersonic  jet,  i.e.,  how  long  is  the  supersonic  region? 

3.  For  (2),  how  do  the  parameters  K,  Ke,  u>  and  y  affect  the  penetra¬ 

tion  length  and  the  relaxation  and  decay  of  the  complex  nonlinear 
wave  interaction  process  to  uniform  conditions? 

4.  Under  what  circumstances  does  a  throat  arise  in  the  jet? 

5.  Are  any  of  the  ideas  associated  w’th  the  one  dimensional  flow 

theory  of  under-  and  over-expanded  plumes  applicable  to  delineate 
the  various  flow  regimes  in  this  two-dimensional  context? 

6.  If  t^e  streamwise  scale  of  the  reflections  is  0(5^),  then  the  ex¬ 
ternal  flow  sees  a  wavy  wall  boundary  condition  if  the  jet  remains 
supercritical .  In  this  case,  the  oasic  slip  line  shape  is  modulated 
by  corrugations  induced  by  wa^e  reflections,  and  the  relaxation 
length  of  the  interaction  process  may  not  be  of  0(6^3).  Are  Poincare 
limit  process  exoansions  of  the  type  given  in  Eqs,  (1)  -  (5) 
adequate  to  treat  the  problem,  or  is  “two-timing"  required? 

During  this  phase  of  the  research,  a  partial  answer  to  some  of  these  questions 
has  been  obtained  in  connection  with  an  integral  mass  conservation  form  of  the 
KG  equation  tor  the  flow  in  the  jet,  which  is 
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l“dy  /  <fi|(Ke»;  +  <1  • 


If  the  jet  slip  line  is  given  by  y  -  s4^  5(x),  where  x  =  6_ly/3x,  then  the 


tangency  condition  is 


$y(x,0)  =  B'(x) 


Noting  that  5(0)  =  0,  (13)  and  (14)  lead  to  the  relation 


5(oo)  =  -  /  (Kg4>-  -  -^2^-  -}  _  dy  . 

-1  X-**> 


Furthermore,  if  we  assert  that  ^  and  ^  +  0  as  x  ►  «,  then  the  slip  line 
pressure  condition 


(x,0)  =  ^-(x,0)  -  w/r  (16) 

implies  if  further,  ambient  pressure  is  achieved  in  the  freestream,  i.e., 

4>-(x ,0)  -►  0  as  x  ->•  oo,  that 

‘P-  *  w/y  as  x  +  oo  (17) 

at  least  on  the  jet  boundary  y  =  0.  If  we  continue  this  in  the  interval 
-1  <  y  <  0,  we  note  that  it  is  a  particular  solution  of  the  K.G.  equation 
satisfying  the  wall  boundary  condition  ^(x,0)  =  0.  Such  a  solution  represents 
an  asymptotic  velocity  distribution  infinitely  far  downstream  which  is  uniform 
across  the  jet  and  is  intuitively  reasonable  on  physical  grounds.  By  contrast 
to  the  limit  x  -oo,  the  second  order  exponentially  small  corrections  for  x  +  <*> 
must  be  determined  from  the  full  K.G.  equation  rather  than  a  Prandtl  Glauert 
approximation.  For  critical  conditions  at  this  location 
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which  is  achieved  by  virtue  of  (17)  when 
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(18) 


Ke  >  ^  (19) 

if  o)  =  1.  Equation  (17)  with  (15)  implies  that 


Tables  2  and  3  as  well  as  Figs.  23  and  24,  summarize  cases  that  can  oc¬ 
cur  by  virtue  of  the  previoi  cormulas  for  K  >  0.  These  cases  are  significant 

for  u$3  application.  Here,  xt  denotes  the  streamwise  position  of  the  throat  or 

minimum  area  section  of  the  jet.  Only  one  configuration  (Case  V)  exists  for 

oj  =  0 .  For  w  =  1  and  u  =  -1 ,  however,  many  configurations  exist,  and  a  diagram 

delineating  the  various  flow  regimes  is  given  in  Fig.  23.  Sketches  of  the  cor¬ 
responding  jet  shapes  are  given  in  Figs.  25  and  26.  Note  that  Case  V  represents 
a  uniform  flow  for  the  jet  which  shears  the  external  freestream,  i.e.,  both 
flows  are  uniform  in  this  case.  Also,  the  indicated  expansion  and  contraction 
of  the  jet  is  consistent  with  one-dinensional  reasoning  for  the  various  cases. 
Note,  however,  that  the  question  of  smooth,  shockless  decelerations  and  tre 
realizability  of  these  inviscid  cases  must  be  resolved  through  experiments  and 
other  analyses  which  could  assess  the  role  of  viscous  entrainment  along  the 
shear  layer.  Computational  studies  could  illuminate  the  question  regarding  the 
existence  of  shocks  in  the  jets. 

In  related  IRSO  work,  the  case  of  Ke,  K  -®  has  been  analyzed, 
involving  decay  processes  associated  with  radiation  of  the  jet  pressure  field  to 
the  external  coflow  through  transmitted  waves.  If  Mg  and  denote  jet  exit  and 
freestream  Mach  numbers,  respectively,  then  relaxation  to  uniform  state  takes 
place  through  a  wave  train  produced  by  n.xi « t i pi e  wave  reflections  from  the  wall 
and  slip  line.  For  the  N*h  downstream  Mach  diamond,  the  analytical  solution  of 
this  problem  indicates  that  the  wall  pressure  Cp  behaves  like 
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TABLE  2  ~  OCCURRENCE 


Case 

No. 

“ 

Jet  Exit 
Supersonic? 

Range  of  1^ 

I 

1 

No 

0  <  *e  <  3r 

la 

1 

No 

Ke-*r 

II 

1 

No 

ill  <  K  < 

Zy  e  Y 

I  la 

1 

No 

K  .  Hi 
e  Y 

III 

1 

No 

Y  « 

IV 

I 

Yes 

—  <  K  <0 
e 

V 

0 

Yes  and  No 

-  «  Ke  <  - 
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OF  CRITICAL  CONDITIONS  AND  THROATS  IN  JET 
(K  >  0) 

(u  »  1 ,0) 


b(-) 

Streamwlse 
Position 
of  Throat 
(x  t) 

Mach 
No. 
at  « 

Monotone 

S1 1  pi  1  ne 
Admissible? 

Remarks 

>0 

<- 

>1 

No 

¥>  («)  *  Y*1 

P  >  P 

e  » 

0 

<- 

>1 

No 

<0 

<» 

>1 

No 

<0 

0* 

-1 

Yes 

<0 

0 • 

<1 

Yes 

>0 

0 

>1 

Yes 

0 

<>1 

Yes 

(•)  *  -Y*1 

p  »  p 

e  " 

'P  *  $  *  0 

Pure  slip  flow 
(No  Perturbations) 

Jet  Exit 
Supersonic? 

Range  of  Kg 

B(-) 

Yes 

<  K  < 
e  y 

i- 

<0 

Yes 

e  y 

<0 

Yes 

<0 

Yes 

v-ir 

■°  i 

i' 

Yes 

! 

-  5^-  <  Ko  <  0 

2y  e 

>0 

Sonic 

>  0 
e 

>0 

NO 

K  >  0 
e 

>0 

Streamwise 

Position 


Mach 
No. 
at  <* 

Monotone 
i  Slipline 
Admissible? 

>1 

Yes 

1 

Yes 

<1 

No 

<1 

NO 

<1 

NO 

<1 

Yes 

<1 

Yes 

Fig.  24  Selection  diagram  for  jet  downstream  criticality  and  throat 
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C  &  -  {1  -  (-1)^  exp  (-N  log  a)}  as  N+« 


where 


M2  a 
e  e 

6  _  REFLECTION 

:  TP  T  "  COEFFICIENT 
e  e 

TP  +  T 


=  /M2"  -  l 
00 


0  =  /TP-TT 
e  e 


LUMPED 

u  _  PRESSURE 

~np  rr  ratio 

l  w*  T  PARAMETER 

00 


If  <P*  is  the  perturbation  potential  in  this  process  related  to  the  exact 
potential  by  the  asymptotic  expansion 


4>  • 

rp  -  X  +  W*{x*,y)  +  0(6Z) 
e 

for  x*  =  x/6,  y  =  (y  -  fib  (x))/6,  M  ,  M  fixed  as  6+0,  then  v*  satisfies  an 

X  c  ® 

initial  boundary  value  problem  for  the  supersonic  Prandtl  Glauert  equation 
giving  the  damped  wave  train  solution 

“■  (x*  +  3ey)  £(-a)°H[x*  +  Bgy  -  (2n+l)eJ 

+  (x*  -  eey)  £(-a)nH  [x*  -  eey  -  (2n+i)pJ  (; 


where 


H(x)  =1  ,  x  >  0 

=  0  ,  x  <  0 


0  <  a  <  1. 


The  problem  of  incompressible  flow  is  also  interesting.  With  a  slight 
change  in  normalizations,  we  have  for  this  case 
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A<£=0  |  -oo  <  X  <  <» 

in  j  -1  <  y  <  0 
=  0  ,  y  >  0 

with  the  boundary  conditions 


^y(X,0)  =  Cj)y  (  X  »  0  ) 

^y(x,0)  =  <fry(x»0) 

^V(x,-D  =  0 

J  mP 


^x(=°.o) 


=  0,  x  <  0 
=  b'(x),  x  >  0 


>Px(x ,0)  -  4>x ( x ,0)  =  c,  x  >  0  , 


where  P  is  the  density,  and  PM  and  U  are  the  ambient  pressure  and  freestream 
velocity  of  the  external  flow.  Cartesian  coordinates  similar  to  Fig.  22  are 
used.  Also,  the  jet  pressure  P.  is  assumed  to  be  slightly  different  from  that  of 
the  freestream  measured  by  the  parameter  5  in  a  simiilar  manner  to  that  utilized 
in  the  previous  developments,  i.e., 

P. 

•pi  =1  +  0)6 

CO 

with  o)  =  -1,  0  or  1,  Furthermore,  the  slip  line  is  given  as  y  =  6b(x),  x  >  0. 
Correspondingly,  we  assume  that  the  jet  velocity  Vj  is  expressed  as 


Finally,  the  external  flow  vector  is  represented  as 


A.  a 

u  ' 


+  5$XK  +  tyyj 


<7(S2) 
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and  the  jet  is 

U  =  (1  +  54>x)^  +  + 

where  t  and  J  are  unit  vectors  in  the  x  and  y  directions,  respectively.  As  in 
the  other  speed  regimes,  two  different  potentials  are  needed  because  of  the 
different  Bernoulli  constants  in  the  jet  and  external  flows.  The  homogeneous 
Neumann  conditions  at  x  =  -»  establishes  a  pressure  level  for  the  jet.  The 
corresponding  inhomogeneous  Neumann  condition  at  x  =  »  is  consistent  with 
equilibration  of  the  pressure  to  the  ambient  level  of  the  external  flow  (P„). 

The  aforementioned  boundary  value  problem  leads  to  the  following 
singular  integral  equation  for  b'(x),  when  Green's  function  representations  for  <f> 
and  'P  are  used  in  their  respective  domains  of  validity  and  the  boundary 
conditions  are  applied: 

/  b'U){tanh  (x-0  -  -Mdi-  =  c(tt  -  i) 

0  x'^  c 

which  can  be  solved  by  Weiner-Hopf conformal  mapping  -  Green's  function 
methods,  or  numerical  procedures.  This  formulation  can  be  used  to  check  the 
wel 1 -posedness  of  the  transonic  problem  and  the  appropriateness  of  the  numerical 
procedures  which  we  will  now  discuss. 

3.1.2  Basic  Ideas  for  the  Computational  Solution  of  the  Problem  of  Fig,  22 

Because  of  the  nonlinearity  of  the  flow  problem  represented  in  Fig.  22, 
numerical  methods  must  be  used.  The  approach  that  has  been  implemented  employed 
a  SLOR  Murman-Cole  type  dependent  difference  scheme  with  Jameson  damping,  in 
which  the  QP  slip  line  pressure  boundary  condition  is  tangentially  integrated  to 
give 


<t>(x,0)  =  v>(x,0)  -  x  +  C  ,  x>0 


(23) 
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where  C  =  ,(>(0,0)  - 
the  relation 


cp(0,0)  is  updated  in  successive  sweeps.  Equation  (23)  and 


'J’y(x.O)  *  *-(x,0)  ,  x  >  0  ,  (24) 

with  the  other  boundary  conditions  indicated  in  the  figure,  define  the  numerical 
problem.  In  this  connection,  a  far  field  in  the  external  field  given  by 

V$  =  0{r~l)  -  0  ,  as  r2  =  x2  +  Ky2  »  (25) 

is  prescribed  on  the  finitely  large  computational  boundaries  P'Q’,  Q'Q  and  QP. 
Compatible  with  (25), 

(j/y  as  x  co 

which  also  is  satisfied  on  the  finite  computational  boundary  PR.  For  ease  of 
calculation  and  programming  logic,  dummy  matrix  entries  have  been  used  in  the 
region  LABP1  in  Fig.  22. 

Details  of  Numerical  Procedure 

In  some  respects,  the  numerical  procedure  is  similar  to  that  discussed 
in  Ref.  5.  Referring  to  the  formulation  in  Fig.  22,  for  interior  points,  the 
transonic  potential  equation  for  the  jet  region  and  external  region  in  diver¬ 
gence  form  is  discretized  by  using  central  differences  when  the  equation  is  ellip¬ 
tic  and  backward  differences  when  it  is  hyperbolic.  Parabolic  and  shock  point 
operators  are  employed  when  required.  Both  the  Neumann  boundary  conditions  and 
the  subsonic  and/or  supersonic  jet  conditions  are  discretized  as  described  in 
Ref.  5.  New  to  this  problen  are  the  interface  conditions  (24)  and  (25). 
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The  former  can  be  written  as 

4>{x  ,0)  _  ^(XjO)  -  v(x  -  ^jjp)  ^jjp  *  (26) 

Here,  v  =  ±1/y  and  cITF  is  a  constant  to  be  determined.  The  quantity  XjTp  is  the 
x-coordinate  of  the  cowl  lip. 

Since  no  mesh  points  are  placed  on  the  line  y=l,  Cjyp  is  determined  by 
extrapolation  of  <j>  from  above  and  $  from  below  the  slip  line  along  x=Xjyp. 

Tnus,  c jjp  =  ^(xitp>1)  -  ^(xitp*1)* 

Equation  (26)  is  inserted  into  the  discretization  of  the  operator  <|>  at 
the  first  '^sh  point  above  the  interface  by  replacing  ^  by 

“  v(xi  "  x 1 TF ^  "  C ITF  * 

Similarly,  below  the  interface  uses  4>i  j+i  replaced  by  <t>-jj+i  +  v(x-j  -  x jyp) 

+  CjTp.  This  type  of  differencing  guarantees  that  Eq.  (26)  is  automatically 
satisfied.  The  relaxation  sweep  is  carried  out  as  in  Ref.  5  from  upstream  to 
downstream,  where  homogeneous  Dirichlet  conditions  are  utilized  in  Station  AB  in 
Fig.  22. 

INITIAL  PARAMETRIC  MALL  JET  STUDIES 

Prior  to  dealing  with  conditions  (23)  and  (25),  the  implications  of  the 
initial  conditions  and  pressure  gradients  in  the  problem  of  Fig.  22  were 
investigated  in  parametric  studies  that  were  conducted  with  an  extended  version 
of  our  original  submerged  wall  jet  code  described  in  previous  sections.  The 
problem  studied  was  in  the  previous  notation: 


T[  ^Ke]  =  [Ke  "  (T+i^xl'Axx  +  *yy  =  0 
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*x(0,y)  3  PHIXO  =  g(y) 

v>(0,y)  =  h(y)  ,  (specified  for  hyperbolic  exits) 

^(x,0)  =  G(x)  ,  x  >  0 

^y(x.-l)  =  0 

with 

v>(x.y)  =  ^pF(y)  as  X  . 

Note  that  in  contrast  to  the  formulation  including  Eqs.  (23)  and  (24),  which 
involves  the  determination  of  the  pressure  on  the  slip  line,  the  formulation 
here  assumes  this  variation. 

Table  4  provides  a  tabulation  for  some  of  the  cases  run  for  this 
problem  where  v>pp  was  assumed  to  vanish  for  these  studies.  Comments  on  the 
qualitative  features  of  these  examples  are  provided  in  this  table. 

From  these  studies,  Figs.  27  and  28  illustrate  typical  results  for 
Ke  =  0.1  for  the  effect  of  the  initial  velocity  g(y)  on  penetration  for  a  submerged 
jet  with  <n-(x,0)  =  0.  Here,  Cp  =  (Cp)/(62^3) ,  and  C*  is  the  critical 
pressure.  Foe  h  =  -1,  the  jet  remains  subcritical.  For  the  h  >  0  cases,  it  is 
evident  that  the  shock  terminating  the  supersonic  region  moves  progressively 
further  downstream  with  supercritically  increasing  exit  velocity. 

Figures  29  and  30  depict  the  effects  of  assigning  a  coflow  pressure 
gradient  along  the  slip  line  ^-( x ,0)  =  1  as  compared  to  the  submerged  jet  case, 
<p-(x,0)  =  0.  I  he  former  is  crudely  analogous  to  (23).  It  is  evident  that  the 
pressure  gradient  has  a  dramatic  influence  in  accelerating  the  flow  and 
prolonging  the  supersonic  zone.  This  fact  could  affect  the  matching  with  the 
Spence  region  for  the  problem  of  Fig.  2ia  in  the  inherent  assumption  that  the 
wave  interaction  zone  there  is  r>(o^). 
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0 
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18 

0 

c 
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19 

0 

0 
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TABLE  4  -  PARAMETRIC  SUMMARY  Of  RUNS 

♦(*•!)  f  ^SUP.  >y(*»0) 

0  0.1  686  0 

0  0.1  118  0 

0  lO'3  764  0 

x  10'3  2652  0 

x  -10*3  2652  0 

0  10-5  0  0 

0  0.1  00 

0  0.1  164?  -2x,(0,l) 

0  0.1  1435  -2x,(0,l) 


-2*. (0.1) 
-2. (1.2) 
0.(2, 3. 8) 


Ct  ,(3.8, 6. 2) 

0. (0,3.8) 

Ci ,(3.8, 6. 2) 


Mostly  elliptic.  Exit  hyperbolic, 
Ter  field  elliptic. 

Completely  elliptic  except  Jet 
exit  which  Is  hyperbolic. 

Jet  exit  supersonic. 

Mostly  hyperbolic  with  elliptic 
fer  field. 

Mostly  hyperbolic  with  elliptic 
far  field. 


Completely  elliptic. 

Larger  supersonic  region  than 
-2,(1 ,2)w1thout  wall. 

0.(2.-) 

Super'onlc  region  smaller  than 
-2,(1. 2)case  with  less  blowlnO. 
0.(2,-) 

Smaller  supersonic  region  than 
case  with  wall . 

*(x,0)  .  C,  In  (3. 8,6. 2). 

Field  similar  to  Case  8  without 
slip  line  except  for  In  (3.8,0). 

Supersonic  everywhere. 

Supersonic  everywhere  (with  wall). 

Subsonic  jet  exit.  Very  small 
supersonic  far  field. 

Supersonic  almost  everywhere. 

Supersonic  almost  everywhere. 

Supersonic  everywhere. 

Subsonic  everywhere. 

Supersonic  everywhere. 

Subsonic  everywhere. 
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Submerged  wall  jet  studies  (continued)  -  initial  condition  effect  on 
extent  of  penetration  of  supersonic  zone,  Kp  =  0.1,  YSONIC  =  position 


I  Rockwell  International 

Scmoc*  C*nt*r 

SC5055.21FR 


O 

■  • 

&  b 


o 

• 

QJ 

L_ 

</> 

• 

to 

» 

0 

0) 

l» 

Q. 

a' 

d 

• 

■r- 

10 

> 

CL 

10 

In 

• 

4- 

▼ 

► 

D 

O 

4-> 

o 

<D 

4— 

,  • 

«- 

4— 

<D 

1  • 

ro 

v> 

-a  o 

t* 

1 

(D  T— 1 

rj 

c  It 

Z>  <D 

C  ^ 

o 

o  * 

r> 

lies  ( 
shape 

J'OlOiNU 


79 


SC5055.21FR 


3.1.3  Results  for  Complete  Problem 

Applying  the  numerical  procedure  described  previously  to  treat  the 
problem  involving  Eqs.  (23)  and  (24),  results  were  obtained  for  illustrative 
cases.  In  this  connection,  a  typical  three-dimensional  relief  plot  of  Mach 
number  as  a  function  of  spatial  coordinates  x  and  y  related  to  the  formulation 
in  Fig.  22  is  shown  in  Fig.  31,  with  y  now  measured  from  the  bottom  wall.  In 
the  legend,  recall  that  the  quantity  ta  is  given  by 


2/3 

where  Pe  is  the  jet  static  pressure,  PM  is  that  of  the  freestream,  and  <j>  '  is  a 
small  parameter  measuring  the  indicated  excess  pressure  ratio.  The  quantity  w 
denotes  a  sign  Index  which  can  be  1,  -1,  or  zero.  In  addition,  the  parameters 
Kfi  and  Kj  are  defined  as 

Ke  =  (1  -  M2)/S2/3  ,  Kj  =  (1  -  M2)/62/3  . 

For  the  case  indicated  in  Fig.  31,  w  =  1,  K  =  1,  and  K.  =  .  In  the  code, 

“  J  Y 

the  jet  exit  potential  <p(-l,y)  was  assumed  zero.  In  contrast  to  the  homogene¬ 
ous  Neumann  conditions  assumed  in  the  derivation  of  Table  1,  this  procedure  has 
been  justified  in  an  earlier  section  of  this  report. 

In  spite  of  the  finite  location  of  the  jet  upstream  boundary,  the 
results  obtained  are  roughly  compatible  with  those  predicted  in  Table  I.  Fig¬ 
ure  31  shows  a  strong  acceleration  of  the  initially  subsonic  jet  flow  embedded 
in  a  subsonic  stream  jet  flow  due  to  choking  associated  with  constriction  of  the 
jet  slip  line  (not  shown).  A  terminating  shock  (supersonic  to  subsonic)  forms-, 
giving  sonic  conditions  Infinitely  far  downstream. 

Figure  32  depicts  a  flow  pattern  associated  with  an  inltally  supersonic 
jet  in  a  subsonic  stream.  The  acceleration  Is  smooth  and  isentropic.  Both 
figures  ypify  wave  Interaction  regimes  that  can  result  in  global  modifications 
of  the  flow  field,  resulting  from  a  change  in  the  Kutta  condition  wake  structure 
and  interaction  with  the  viscous  wall  jet  on  upper  surface  blown  airfoils. 
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3 . 2  Solution  of  Blown  Airfoil  Problem  and  Relevant  Flow  Structures 

In  this  section,  in  contrast  to  the  detailed  jet  studies  previously 
given,  we  change  focus  and  concentrate  on  various  formulations  and  analyses 
which  are  intended  as  a  critical  assessment  of  different  aspects  of  the  jet  flap 
formulation  given  in  Ref.  1  and  its  extension  to  USB  for  incompressible  and 
compressible  flows.  Whereas  the  previous  discussions  have  concentrated  on  wall 
jet  flows  relevant  to  the  upper  surface  blown  problem,  this  discussion  will 
treat  other  formulational  issues.  In  particular,  as  a  basis  for  the  thin 
airfoil  and  small  disturbance  models,  the  thin  jet  approximation  is  described  in 
Section  3.2.1  in  terms  of  a  systematic  asymptotic  expansion  procedure  for 
incompressible  flows.  This  study  is  used  as  a  basis  for  obtaining  boundary 
conditions  across  the  jet.  The  generalization  of  these  developments  to 
compressible  flow  is  straightforward  and  therefore  not  provided.  Section  3.2.2 
provides  a  discussion  of  the  trailing  edge  region  from  the  viewpoint  of  a 
nonuniformity  of  the  thin  jet  theory  as  well  as  the  generalized  Kutta  condition 
for  USB.  Inherent  in  this  aspect  is  the  geometry  of  the  dividing  streamline  at 
the  trailing  edge,  which  is  a  necessary  condition  for  the  determination  of  the 
jet  sheet  free  boundary.  Arguments  are  provided  to  substantiate  tangency  to  the 
upper  surface,  providing  the  jet  stagnation  pressure  is  greater  than  that  of  the 
external  flow.  Also  indicated  is  how  the  small  disturbance  jet  flap  formulation 
of  Ref.  5  is  modified  with  USB.  Finally,  Section  3.2.3  gives  resuUs  from  a 
computational  solution  based  on  the  inviscid  small  disturbance  formulation 
indicated  in  Fig.  20b.  In  this  section,  transonic  USB  airfoils  are  analyzed  and 
comparisons  are  made  with  experiment.  Factors  associated  with  the  discrepancies 
are  considered  and  refinements  are  proposed  to  improve  the  realism  of  the  model, 
leading  to  the  discussion  of  the  viscous  model  given  in  Section  4. 


3.2.1  Thin  Jet  Theory 

As  an  essential  ingredient  of  a  small  disturbance  formulation,  the  jet 
structure  Is  developed  In  this  section  for  purposes  of  specification  of  the 
boundary  conditions.  In  particular,  it  will  be  shown  how  the  Spence  theory  of 
Ref.  1  can  be  derived  and  refined  from  a  systematic  approximation  procedure. 
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Referring  to  Fig.  33,  a  section,  of  the  jet  is  shown.  A  curvilinear 
coordinate  system  is  embedded  in  the  jet  as  indicated.  The  lines  n  =  constant 
are  parallel  to  a  reference  line  (the  g  axis)  which  is  the  (flux)  center  line  if 
the  jet  is  symmetric,  or  the  wall  for  a  wall  jet  case,  in  this  coordinate 
system,  the  lines  £  =  constant  are  normal  to  £  axis.  In  what  follows,  the 
incompressible  case  will  be  discussed.  The  generalizations  to  compressible  flow 
are  straightforward. 


Within  the  indicated  coordinate  system,  the  exact  equations  of  motion 


are 


Continuity 


—  q  +  -- 


h  q  =  0 
en  n 


g  -  Momentum 


n  sg 


qcq 


q  H£  +  XDi»i-4 

0  3n  h  ph 


n  -  Momentum 


n  3g  %  8n  n  3n 


1  3P 
p  3n 


(27a) 


(27b) 


(27c) 


where  h,  the  metric  coefficient  when  related  to  the  differential  arc  length  in 
Cartesian  (x,y)  coordinates  is 

dx2  +  dy2  =  h2(g,n)dg2  +  dn2  {28) 

b(5»n)  =  1  -  n/R(c) 

with  R(g)  being  the  radius  of  curvature  which  is  shown  positive  in  Fig.  33. 

To  obtain  an  approximate  incompressible  set  of  equations  prototypic  of 
the  compressible  case,  the  thin  jet  limit  is  considered.  The  characteristic  jet 
thickness  Is  shown  in  Fig.  33,  where  the  jet  boundary  is  denoted  as  n  *  tb(g). 
The  quantity  t  thus  represents  a  characteristic  slope  of  the  jet  free  boundaries 
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Fig.  33  Section  of  jet  and  curvilinear  coordinate  system. 
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which  is  of  the  same  order  as  the  jet  thickness  in  units  of  its  radius  of 
curvature. 

We  now  define  a  thin  jet  limit 


t  *  0,  £.n*  =  h/t  fixed 


(29) 


where  the  boundary  layer  coordinate  n*  is  introducted  to  keep  the  jet  slip  lines 
in  view  in  the  limit  process.  In  (29),  the  appropriate  representations  to  yield 
a  nontrivial  structure  are 


9r(5»n;t)  , 

-  -g - *  ft  Ug(C.n*)  +  /T" UjU.n*)  +  ••• 


(30a) 


>=/7v0  +  t3/V 


(30b) 


pH 


2  =  P0  +  TP1  + 


(30c) 


where  U  is  some  typical  freestream  velocity.  The  orders  were  selected  to  give 
the  "richest"  possible  set  of  equations  and,  consistent  with  this,  produce 
forcing  terms  in  the  equations  for  the  second  order  quantities.  These  orders 
are  consistent  with  the  massless  momentum  source  mode  of  Spence\ 


Substitution  of  (30)  in  (27)  and  equating  like  orders  gives  the 
following  equations  for  the  approximate  quantities: 


fV  n 

35  3n* 


(31a) 


auj. 

Jo  W 


+  V 


3UC 

0  87? 


(31b) 


u0 

m 


3P0 

3n* 


(31c) 
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3U1 

w + 

tl 

<ro  |<ro 

1  3 

3u. 

1  +  V 

3U. 

A  X 

3U0 

X.  _ \J 

uov 

uo  W 

-+  vo 

3^  + 

35  1 

3n*  1 

"  R 

2uoui 

3Pi 

=  -  “ 

3Vq 

O 

> 

<T> 

u0n* 

R 

+  3n* 

0  35  ‘ 

V0  3?‘ 

'  -RT- 

3Uq  3P0 


TT  u0  35 


H 


(32a) 


(32b) 


(32c) 


The  appropriate  boundary  conditions  involve  statements  regarding  the 
fact  that  the  jet  boundaries  are  streamlines  and  that  the  static  pressure  is 
continuous  across  the  slip  lines.  The  upper  and  lower  slip  lines  Su  and  -  0 
are  given  by 

su  =  n  -  -rb0U)  -  AjU)  =  0 

=  n  +  xb^U)  +  T2b^(5)  =  0 

where  symmetry  has  been  assumed  to  leading  order. 

Based  oi.  the  foregoing  discussion,  the  condition  that  jet  boundaries  S 
are  streamlines  is 


<$  •  VS  =  0 

where  $  =  (q  .q^).  Substitution  of  the  expansions  (30)  into  this  relation  gives: 

V5'b0}  =  b0U) 

vQU.-b0)  =  -bbU)  u0U.-b0) 


87 

C347 1 A/ j  bs 


SC 5055. 21 FR 


vn(?,bn)bn  aun 

v!^«bo)  =  r  +  bibo  aTr*-  ^,bo^  (5*bo)  +  ui^,bo^bo 

3  V 

-  (e.b0)  *  blU()(?,b0) 


(33c) 


vf)(?  »"bn)bn  3urj 

vl(c.-b0)  =  r  +  blb0  Tn*  ^’“b(P  '  ul^‘-b0)b0 


av, 


+  bl  an*  (5»-b0)  “  Vo^'V 


(33d) 


The  other  conditions  Involving  continuity  of  pressure  are: 

p(U)(e) 

P0(«*bo) 

Pl(c*b0)  =  -  h^c)  -g^5F  (?»b0) 


(34a) 

(34b) 


where  we  assume  for  the  present  argument  that  the  external  pressure  field 
is  prescribed.  Similar  pressure  conditions  hold  for  the  lower  slip  line 
boundary.  If  6  is  a  characteristic  flow  deflection  angle  of  the  order  to  the 
airfoil  thickness  or  angle  of  attack,  then  T  «  6  has  also  been  implicitly 
assumed. 


Solutions 

First  Order  Theory 

We  introduce  the  zerotb  order  stream  function  given  by 


14  - 

n  " 


-  v. 


u 


0 


and  employ  the  following  transformation  for  the  Independent  variables 
U.n*)  *  U»*)  • 


(35) 
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Under  (35),  the  differential  operators  map  as  follows: 


J_  ,  i_  _  y  — 

35  35  0  34> 


3n*  =  U0  3* 


and  Eqs.  (31)  become 


3U0  3U0  3V0 

8?  0  dtp  U0  3^  U 

(36a) 

(36b) 

^0  =  _  3£ 

R  3ij) 

(36c) 

The  general  solution  for  Eqs.  (36)  is 


uo  =  u0(») 

(37a) 

fo 5  k<«>  -sxzrl  uo<*'>  <*** 

(37b) 

v0(5 •'i')  =  b‘(5)  u0U)  , 

(37c) 

k(5)  is  the  pressure  on  the  centerline.  The  reference  (centerline)  condition 
(v0(5>0)  =  0)  or  the  conditions  at  the  edges  of  the  jet  (Eq.  (33b))  show  that 
bg  =  0  and  bQ  =  1,  if  bQ  (0)  =  1.  Thus,  in  terms  of  an  arbitrary  initial 
parallel  flow  profile  Ug  =  f(n*)  at  5  =  0,  we  have 

n* 

i>  =  J  f(n)  dn 

0 

(38a) 

Therefore,  ug  =  f(n*(<>))»  and 

(38b) 
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P0U.*)  =  MO  -  ]j||y  A  (n*(0)  d* 

n*  1 

=  MO  -  *  /  f2dn*  =  qu(?)  +  tj  /  f2dn*  (38c) 

0  n* 

(MO  =  qu  +  ^  /  f2dn*)  • 

n* 


Discussion 

Equations  (38)  describe  a  parallel  flow  jet.  The  total  jump  in 
pressure  across  the  jet  from  (38c)  is 


[P03  =  PgU.O  -  Pq(c.-I)  =  -  |  Jlf2(n*)  <V  (39) 

which  agrees  with  the  Spence  model.  It  should  be  noted  that  in  contrast  to  the 
latter,  no  assumption  regarding  irrotationality  is  required  to  obtain  (39), 
which  is  also  in  contrast  to  the  results  of  previous  workers.  The  radius  of 
curvature  of  the  jet  is  approximately  R  upstream  of  the  trailing  edge,  which  in 
turn  is  approximately  given  by  that  of  the  blown  upper  surface.  Downstream  of 
the  trailing  edge,  R  is  determined  from  applying  (39)  to  the  determination  of 
the  flow  outside  of  the  jet.  Upstream  of  the  trailing  edge,  the  wall  pressure 
is  determined  by  (39),  since  R  is  known  and  is  given  by 


PU.-l)  =  quU) 


H 


f2(n*) 


dn* 


(40) 


In  Ref.  6,  the  appropriate  second  order  theory  corresponding  to  the 
expansions  (30)  are  derived.  Also  indicated  are  various  nonuniformities  that 
can  occur,  such  as  those  near  the  trailing  edge.  These  reveal  the  need  for  other 
expansions  in  the  neighborhood  of  the  nonuniformities. 

Also  in  Ref.  6,  a  small  deflection  specialization  of  the  previous  thin 
jet  theory  is  considered.  Therein,  a  systematic  patching  procedure  is  utilized, 
employing  a  "blending  layer"  at  a  vertical  distance  of  the  order  of  the  slip 
line  and  airfoil  deflection  angle  6  as  6  +  0,  The  blending  layer  has  been  used 


90 

C3471A/J  bs 


, '  „ 


Rockwell  International 

Science  Center 


SC5055.21FR 


to  validate  the  usual  Taylor's  series  transfer  of  boundary  conditions  employed 

to  define  the  outer  flow.  A  treatment  of  similar  blending  layers  is  discussed 
23 

in  Cole  in  connection  with  incompressible  flows  around  unblown  bodies  of 
revol ution. 


3.2.2  Trailing  Edge  Behavior 
Incompressible  Flows 

Defining  the  complex  potential  function  as 

F(z)  =  $(x,y)  +  H(x,y) 

where  (x,y)  is  the  local  coordinate  system  shown  in  Fig.  34  and  z  =  x+iy,  the 
local  "corner  flow"  solution  to  within  a  dimensional  multiplicative  constant  in 
the  lower  external  region  A' OB  where  B  and  v  are  the  dividing  streamline  and 
trailing  edge  angles  shown  is  given  by 

P  =  e-iir(ir+v)/B  _,tt/B 

which  implies  that  the  square  of  the  resultant  velocity  q  is 

|F'|2  =  q2^r2^  ^  .  (41 ' ) 


The  initial  conditions  for  determination  of  the  jet,  which  follow  from  the 
requirement  that  the  pressure  and  flow  angle  be  continuous  across  the  dividing 
streamline  near  the  trailing  edge,  were  studied  for  incompressible  flow. 

Referring  to  Fig.  34,  with  the  dividing  streamline  denoted  as  OB,  we  signify  the 
trailing  edge  angle  as  v  and  the  angle  that  OB  makes  with  the  upper  surface  AO  as  u> 
and  that  with  the  lower  surface  as  b.  The  details  of  this  discussion  are  given 
in  Ref.  6  which  shows  that  a  =  7t  and  b  <  n  is  the  only  viable  possibility  if  the 
stagnation  pressure  Pg  above  AOB  is  greater  than  that  below  it,  which  is 
signified  by  Pg,  where  +  and  -  superscripts  signify  conditions  above  and  below 
the  slip  line  AOB,  respectively,  and  0  refers  to  stagnation  conditions. 

Denoting  u  as  the  flow  speed  along  the  slip  line,  Bernoulli  implies 
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Here  also,  u"  =  0.  These  results  quantify  the  slip. 

If,  on  the  other  hand,  P*  <  P‘,  then  w  <  tt  and  e  =  n  is  the  only 
possibility  with 


for  a  hypothetical  case  in  which  P‘  «  P*  in  a  real  physical  flow  (fluctuating 
above  and  below  equality),  then  a  tri-stable  configuration  could  evolve  which 
would  oscillate  between  configurations  in  which  the  dividing  streamline  is  tan¬ 
gent  to  the  upper  and  lower  surface  or  bisects  the  trailing  edge  angle.  Similar 
arguments  have  been  applied  to  treat  the  conditions  of  the  stream  sheet  at  the 
trailing  edge  of  unblown  incompressible  three-dimensional  wings  in  Ref.  24. 

Nonuniformities  of  Second  Order  Approximations 

From  (41),  we  find  that  another  nonuniformity  associated  with  second 
order  solutions  discussed  in  Ref,  6  is  related  to  the  following  behavior 

q^  %  f( 2tt/b)-3  a$  p  +  0 

where,  we  distinguish  the  following  possibilities  as  ^ 

0)  +  00  »  B  >  or  v  < 

( i  i )  (l|j+0,&<-|^Qrv>-j 

(iii)  q^  finite  *0,3=|£orv=-| 

Case  (i)  is  the  most  practical  situation  and  will  necessitate  an  inner  solution 
for  the  transition  layer  to  join  the  wall  and  free  jet  flows.  This  aspect  is 
discussed  in  Ref.  6.  For  the  jet  structure  considered  here,  we  briefly 
investigate  the  vorticity  which  can  be  shown  to  be  given  by 


93 

C347 lA/j  bs 


SC5055.21FR 


-3/2 

(U  =  T 


4,  -1/2  , 

a>0  +  T  uij  + 


where  as  an  Illustration,  if  ug  =  C,  then. 


(Og 

“1 


0 


3U1  uo  c 

w  '  f  =  ■  R[uy  =  constant  *  0  • 


Because  of  the  constant  initial  velocity  profile,  a  non-zero  vorticity  is  intro¬ 
duced  by  the  body  curvature.  By  contrast,  a  potential  vortex  over  a  circular 
cylindrical  surface  would  have  had  a  linear  Initial  profile  to  produce  an 
irrotational  flow.  Note  for  arbitrary  Initial  profiles,  uQ(0,o*)  =  f(n*)  and 
thus,  u)q  =  f'(n*)  *  0. 


Compressible  Trailing  Edges 

Consider  again  the  configuration  of  Fig.  34.  Here  we  analyze  first  the 
case  where 


PS>  P0 


With  the  usual 


isentropic  r< nation 


P 

P 


+ 

V 

0 


(P+  =  static  pressure  on  upper  side) 


the  local  Mach  number,  M,  can  adjust  in  a  continuous  way,  so  that  the  flow 
recompresses  smoothly  from  B  to  0.  At  0,  M  is  single  valued  =  Mj  and  adjusts 
itself  such  that  the  static  pressure  P+  equals  Pg,  In  accord  with  the  isentropic 
relation,  assuming  that  the  flow  stagnates  at  0  on  the  lower  side.  The  only  way 
this  can  be  realized  is  with  the  u  >  tt  and  g  <  n  arrangement.  In  some  respects 
the  configuration  resembles  flow  over  a  solid  wall  expansion  corner.  This  gives 
rise  to  an  expansion  fan  interacting  with  a  sonic  line  from  the  corner.  If  the 
solid  configuration  were  representative  of  this  flow  with  the  free  boundary  slip 
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line,  compression  waves  would  reflect  off  the  sonic  line  and  form  a  shock 
envelope  which  would  be  necessary  to  recompress  the  flow  from  an  overexpanded 
value  be’ow  the  critical  to  the  PZ  level.  Additional  reflections  can  be 
produced  from  the  upper  slip  line  CD.  This  discontinuous  transition  leads  to  a 
multivalued  pressure  at  0.  The  continuous  and  discontinuous  processes  are 
illustrated  schematically  in  Figs.  35b  and  35c,  respectively.  In  Fig.  35b,  the 
recompression  takes  place  on  the  line  30.  In  Fig.  35c  it  occurs  on  00'0"0'" 
signifying  the  confluence  of  multiple  states  at  0.  Here,  the  dashed  line 
element  0"0'"  signifies  a  shock  jump.  Experimental  data  strongly  suggests  that 
-Vj  in  the  incompressible  case,  the  configuration  with  u  =  it  and  0  <  it  is  the 
mo:'  probable  situation.  Presumably,  a  more  rigorous  argument  to  support  this 
conjecture  would  rely  on  some  sort  of  stability  analysis.  Specific  criteria  to 
determine  whether  the  discontinuous  or  continuous  configuration  occurs  in  a 
given  case  is  an  open  question  that  could  benefit  from  further  study. 

For  Pg  =  P“,  t.ha  dividing  streamline  would  again  bisect  the  trailing 
edge,  since  the  flow  in  the  immediate  vicinity  would  be  incompressible  and  the 
reasoning  in  the  previous  section  would  apply.  For  the  improbable  case  of 
Pg  <  Pg,  the  configuration  with  u>  <  n  and  0  =  tt  wculd  be  applicable  as  for  the 
incompressible  situation. 

External  Flow 

For  distances  large  compared  to  the  jet  width,  the  fine  structure  of 
the  jet  is  important  only  insofar  as  it  provides  matching  conditions  to  the 
irrotational  "outer"  flow  field.  In  incompressible  flow,  this  external  outer 
flow  can  be  determined  by  thin  airfoil  theory.  At  transonic  speeds,  small 
disturoance  theory  is  appropriate  for  this  region.  Details  of  the  asymptotic 
matching  procedure  have  been  discussed  for  Incompressible  flow  in  Ref.  6.  Based 
on  these  developments  and  the  earlier  ones  for  thin  jets  in  this  paper,  the 
boundary  conditions  for  the  outer  flow  in  the  incompressible  and  transonic  c.;-t.s 
for  the  jet  flap  and  upper  surface  blowing  are  now  .ndicated. 
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Referring  to  Fig.  36,  the  equation  of  an  airfoil  can  be  given  as 


y  =  5f(x) 
and  the  jet  is 

y  =  69(x) 

where  6  is  the  thickness  ratio  of  the  airfoil,  and  f  is  the  upper  or  lower  surface 
and  involves  the  angle  of  attack  which  is  assumed  to  be  of  the  same  order  of 
6.  Considering  a  small  disturbance  approximation,  we  obtain 


d2y 

Hx* 


[i  *  (r  ]  11+625 


^-,7737?  !  «”<*>  • 


Letting  the  "outer"  expansion  pressure  coefficient  be  represented  as 


P  -  P 


=  Sp(x,y)  +  ••• 


then  by  virtue  of  a  generalization  of  (39) 


Cp(x,0)]  =  -  C0-g"(x)  =  -  2[4>x] 


where 


cj  5  p/  =  0(1) 


(43a) 


and  <j>  is  a  perturbation  potential. 


Equation  (42)  is  the  relation  used  in  conjunction  with  the  jet  tangency 
boundary  condition 
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4>y(x»0)  =  9 ‘  (x)  (43b) 

and  the  airfoil  boundary  conditions  to  determine  the  external  flow  field.  These 
relations  coincide  with  the  relations  derived  by  Spence.  They  can  be  generalized  for 
transonic  flow  by  placing  the  p  inside  the  integrand  in  (43a). 

Upper  Surface  Blowing 

To  treat  conditions  on  the  blown  part  of  the  airfoil,  Eq.  (42)  can  be 
applied  by  approximating  the  radius  R  by  (f")-1  to  obtain  the  wall  pressures, 
and  by  using  the  airfoil  and  jet  boundary  conditions  to  determine  the  upper  slip 
line  jet  pressures. 

From  the  thin  jet  theory  derived  earlier,  it  can  be  seen  that  rotation¬ 
al  flow  produces  the  same  pressure  jumps  across  the  jet  in  the  dominant  approxi¬ 
mation  as  the  irrotational  Spence  models.  Correspondingly,  It  can  be  shown  that 
to  within  factors  involving  the  density,  qualitatively  similar  results  are  ob¬ 
tained  for  transonic  flow.  Another  important  aspect  of  the  asymptotic  represen¬ 
tations  derived  here  is  that  they  lead  to  higher  approximations  for  the  struc¬ 
ture  of  the  jet  and  external  flow  which  can  be  systematically  obtained. 

Finally,  the  analytical  solutions  described  above  in  Ref.  6  allow  the  systematic 
assessment  of  the  effects  of  initial  vorticity  and  skewness  which  are  Inacces¬ 
sible  to  other  theories. 

3.2.3  Results  and  Discussion  for  Transonic  Upper  Surface  Blowing 

A  successive  line  overrelaxation  (SLOR)  scheme  within  a  Karman  Guderley 
framework  has  been  used  to  compute  the  flow  field  over  an  upper  surface  blown 
airfoil.  On  the  blown  portion,  the  jump  conditions  across  the  jet  are  deter¬ 
mined  by  the  asymptotic  results  given  in  previous  sections,  i.e.,  Eqs.  (42)  and 
(43b).  Providing  that  the  region  is  not  too  close  to  the  jet  exit  or  trailing 
edge,  the  streamwise  gradients  can  be  neglected  in  the  entropy  and  velocity  com¬ 
ponent  parallel  to  the  wall.  Away  from  these  regions,  the  pressure  gradient 
perpendicular  to  the  streamlines  <s  balanced  by  centrifugal  force.  For  the 
region  near  the  jet  exit,  these  assumptions  become  invalid.  Here,  the  scale  of 
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the  gradients  in  the  streamwise  direction  becomes  important,  principally  due  to 
the  influence  of  wave  interactions  with  the  slip  line.  Similar  fine  structures 
occur  near  the  trailing  edge  where  the  flow  can  stagnate  on  the  unblown  side, 
depending  on  the  ratio  of  the  stagnation  values  above  and  below  the  dividing 
streamline.  For  incompressible  flow,  the  previous  sections  have  discussed  the 
tri -stable  equilibrium  at  the  trailing  edge  corresponding  to  the  value  of  the 
stagnation  pressure  ratio,  which  leads  to  the  dividing  streamline  leaving  tan¬ 
gent  to  the  upper  surface  if  this  is  greater  than  unity.  Consistent  with  the 
previous  discussion,  the  appropriate  generalization  to  transonic  flow  was  as¬ 
sumed  also  to  be  this  arrangement  for  a  single  valued  pressure  without  a  shock 
in  that  location.  This  assumption  has  been  altered  to  assess  the  sensitivity  of 
the  flow  to  the  dividing  streamline  angle.  In  this  connection,  surface  pres¬ 
sures  for  the  dividing  streamline  bisecting  the  trailing  edge  angle  (as  it  would 
in  incompressible  flow)  were  compared  with  those  for  the  tangent  arrangement. 
Based  on  these  studies,  significant  differences  are  anticipated  only  for  large 
incidences  and  trailing  edge  angles. 

Typical  results  obtained  from  the  computational  model  are  shown  in 

Fig.  37  in  which  the  flow  over  a  thick  airfoil  designed  at  Rockwell's  Columbus 

Aircraft  Division  (CAD)  was  analyzed  with  the  SLOR  code.  Here,  the  pressures 

for  various  values  of  the  blowing  coefficient,  Cj,  are  compared  with  those 

for  the  unblown  case  at  a  freestream  Mach  number  M  =  0.703  and  anale  of  attack 

00  w 

a  =  0°.  Substantial  lift  augmentation  is  evident  for  blowing.  Also  evident  is 
the  associated  rearward  motion  of  the  shock  with  Increased  blowing  and  sectional 
loading  as  if  the  Incidence  is  increased. 

Further  parametric  studies  are  provided  in  Fig.  38  which  indicate  the 
effect  of  parallel  displacement  of  the  slot  Xj  (in  units  of  the  chord),  on  the 
chordwise  pressures.  Three  positions  of  the  slot  Xj  =  0.5,  0.65,  and  0.8  are 
shown.  No  systematic  trend  in  the  bl own  pressures  Is  exhibited  on  this  airfoil 
with  downstream  slot  movement  for  fixed  Cj.  Evident,  however,  is  a  slight  in¬ 
tensification  of  the  terminating  shock  with  slot  downstream  motion  although  its 
position  remains  unaltered.  Despite  the  limitations  of  the  model  to  describe 
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the  fine  structure  of  the  jet  exit  region,  a  small  suction  peak  which  has  been 
observed  in  experiments  is  exhibited  in  this  vicinity  for  Xj  =  0.5.  In  Fig.  39, 
the  corresponding  increase  in  lift  coefficient  CL  with  slot  downstream  movement 
is  also  shown,  as  is  the  increase  in  the  size  of  the  supersonic  region. 

In  Fig.  40,  the  increase  of  lift  with  blowing  coefficient  as  well  as 
size  of  the  supersonic  region  is  quantified. 

Tests  of  the  adequacy  of  the  foregoing  model  to  simulate  realistic  tran¬ 
sonic  USB  airfoil  flows  have  been  inhibited  by  the  lack  of  suitable  experimental 
data.  Information  exists  only  for  highly  three-dimensional  configurations, 
large  thickness,  or  incidence  in  ranges  beyond  the  validity  of  the  assumptions 
of  small  disturbance  theory.  Another  restriction  is  the  unavailability  of  the 

associated  geometric  data  and  flow  diagnostics  accompanying  the  tests.  The 

pc 

results  of  Yoshihara  and  his  coworkers  were  useful  in  this  connection  and 
allowed  us  to  compare  the  jet  flap  specialized  version  of  the  USB  theory  with 
the  data  in  Ref.  21.  For  the  simulations  described  in  this  paper,  tests  per¬ 
formed  by  N.C.  Freeman  at  NPL  on  a  USB  modified  6%  thick  RAE  102  airfoil,  and 
described  in  Ref.  26,  appear  to  have  the  most  suitable  results  for  comparison  at 
present.  Unfortunately,  the  angle  of  attack  associated  with  the  NPL  data  is  6°, 
which  is  marginal  for  the  application  of  a  small  disturbance  model. 

Figure  41  indicates  comparisons  of  chordwise  pressures  for  various 
values  of  C.;.  Also  shown  are  schlierens,  indicating  the  associated  flow  field 
structure.  Turning  to  the  Cj=0  results  (Part  (a)),  massive  shock  induced  sepa¬ 
ration  is  indicated  and  is  apparently  initiated  at  the  downstream  limb  of  the 
lambda  shock  on  the  upper  surface.  This  is  reflected  in  the  classical  erosion 
of  the  suction  plateau  and  is  responsible  for  the  indicated  disagreement  between 
the  inviscid  computational  results  and  the  data.  For  these  tests,  nominal  tan¬ 
gential  blowing  with  a  slot  height  of  0.07 %  of  the  chord  was  used.  The  slot  loca¬ 
tion  Is  15%  downstream  of  the  nose.  The  Mach  number  M  immediately  above  the  slip 
line  at  the  slot  is  approximately  1.29  for  both  Cj's  indicated.  For  Cj  =0.017 ,  the 
slot  Mach  number  Me  has  been  estimated  as  1.79  and  for  Cj  =0.048,  Mg  «  2.36. 
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Comparison  between  theory  and  ixperiment  in  Part  (b)  of  Fig.  41  indicates 
reduced  discrepancies  on  the  upper  surfc.ee  associated  with  the  limited  separa¬ 
tion.  In  Part  (c),  the  agreement  is  correspondingly  further  improved. 

To  achieve  adequate  realism,  it  is  important  to  discuss  factors  respon¬ 
sible  for  the  disagreements.  One  feature  not  captured  by  the  USB  simulation  is 
the  pressure  spike  at  the  slot  location.  Based  on  the  slot  size,  the  streamwise 
scale  for  this  phenomenon  is  at  least  an  order  of  magnitude  greater  than  the 
characteristic  wavelength  of  a  Mach  diamond  pattern  in  the  wall  jet.  These 
fluctuations  may  not  be  resolvable  with  conventional  pressure  tap  arrangements  for 
the  thin  slot  employed  in  the  tests.  If  a  rough  model  of  a  coflowing  inviscid 
supersonic  wall  jet  over  a  flat  plate  is  used  to  describe  the  flow  near  the  slot, 
the  approach  to  a  final  steady  state  may  be  damped  oscillatory  or  monotone,  depend 
ing  on  whether  the  reflection  coefficient  R  which  is  given  by 

R  =  $  (44) 

where 


A  =  M2fj/M2(5£ 


B  =  /RTT-T 


Bq  +  /wprr 

e  e 


is  respectively  positive  or  negative. 


The  relaxation  length  L  to  achieve  the  downstream  pressure  in  units  of 
the  exit  height  is  of  the  order  of  £n  R,  which  can  be  approximately  5  to  50  in  the 
present  case,  depending  on  the  accuracy  of  the  estimate  for  MQ.  Note  in  this 


connection  that 


R  <  0  for  1  <  M  <  —  ,  and  Me  M  <  » 


M2 


R  >  0  for  <  M  <  M 
8e  e 
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For  the  submerged  case,  R  ->■  1,  (Me  »  M),  and  the  Prandtl  periodic  pattern  is 
obtained,  with  no  radiation  of  energy  to  the  external  flow. 

These  facts  suggest  that  one  factor  that  may  be  responsible  for  the 
observed  spike  is  the  internal  decay  process  in  the  jet.  If  transonic  effects  and 
wall  curvature  are  accounted  for,  the  presence  of  "ballooning"  and  throats  in  the 
jet  may  also  be  contributory.  We  have  discussed  such  phenomena  in  connection  with 
submerged  transonic  wall  jets  in  Ref.  5  and  earlier  in  this  report  and  have 
reported  analogous  results  for  the  coflowing  case  in  previous  sections.  Selection 
rules  in  terms  of  Me  and  M  for  the  existence  of  throats  in  the  jet  near  field  have 
also  been  given  previously  and  are  based  on  an  integral  form  of  the  Karman 
Guderley  eqiation.  A  rough  sketch  of  the  wave  system  that  could  explain  the 
spikes  in  Figs.  41b  and  41c  is  shown  in  Fig.  42.  Yet  another  phenomenon  that  would 
have  a  similar  wave  pattern  would  be  a  slight  upward  motion  of  the  jet  due  to 
viscous  mixing  or  a  misalignment  with  the  surface  tangent  at  point  A. 

Turning  now  to  the  discrepancy  of  the  values  shown  on  the  rear  surface 
(downstream  of  0.5c)  in  Fig.  41c,  we  note  that  in  spite  of  the  obvious  elimination 
of  separation,  a  thick  viscous  wall  jet  is  present.  Downstream  diffusion  will 
affect  the  appl i cation  of  the  Spence  relation  on  the  blown  portion  as  well  as  the 
shock  jump.  In  view  of  the  wall  jet  thickness  shown  on  the  schlierens,  this 
factor  appears  to  oo  more  significant  than  shock  obliqueness  at  its  foot.  A  near 
term  refinement  being  implemented  employs  a  coupled  inviscid-viscous  model  which 
uses  second  order  boundary  layer  corrections  to  the  Spence  boundary  conditions 
accounting  for  axial  gradients  of  the  displacement  and  momentum  thickness.  Once 
these  refinements  are  incorporated,  systematic  optimization  among  separation 
suppression,  wave  drag  minimization,  and  supercirculation  control  will  be 
possible.  It  is  envisioned  that  the  design  techniques  contained  in  Refs.  27-29 
will  augment  this  capability  by  providing  methods  to  modulate  shock  formation  in 
concert  with  the  blowing  effects. 
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Fig.  42  Inviscid  wave  pattern  associated  with  slip  line  ballooning  and  jet 
throat  formation. 
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4.0  VISCOUS  EFFECTS 

4.1  Preliminary  Remarks 

A  schematic  of  the  blown  airfoil  is  shown  in  Fig.  13,  where  the  shaded 
region  indicates  the  boundary  layer  and  wake,  P  is  the  stagnation  point,  and  S  is 
the  slot  location.  In  the  research  performed,  a  computational  model  has  been  de¬ 
veloped  to  treat  the  wall  jet,  wakes,  and  boundary  layers,  which  without  loss  of 
generality  have  been  assumed  laminar.  For  treatment  of  the  latter,  integral 
methods  have  been  discarded  due  to  their  highly  empirical  nature,  lack  of  con¬ 
sistency,  and  questionable  applicability  to  blowing  cases.  In  fact,  despite  the 
unavailability  of  general  accuracy  criteria,  it  is  well  known  that  even  for  un¬ 
blown  incompressible  decelerating  laminar  flows  in  adverse  pressure  gradients, 

30-32 

integral  methods  such  as  those  due  to  Pohl hausen  and  others  produce  results 
that  are  increasingly  inaccurate  as  the  separation  point  is  approached.  For  pur¬ 
poses  of  the  present  tangential  injection  study,  this  limitation  is  unacceptable. 

On  the  other  hand,  full  elliptic  Navier-Stokes  solvers  such  as  those  employed  by 
Diewert  have  not  been  selected  at  this  time  because  of  their  lengthy  computational 
nature  and  associated  need  to  resolve  grid  generation  issues  peculiar  to  the 
present  blowing  problem.  By  contrast,  we  believe  that  the  most  appropriate  model 
is  a  zonal  one,  in  which  the  inviscid  region  is  treated  by  an  iterative  scheme  such 
as  Murman-Cole  SLOR  or  any  of  the  more  modern  techniques,  and  the  viscous  regions 
are  treated  by  some  parabolic  solver  such  as  Keller's  box  scheme.34  A  somewhat  more 
sophisticated  approach  involving  a  parabolized  or  thin  layer  Navier-Stokes  algorithm 
represents  a  future  direction  for  refinement.  The  utilization  of  a  finite  differ¬ 
ence  approach  for  the  viscous  layers  avoids  the  need  for  assumption  of  velocity 
profiles  as  in  the  integral  methods.  This  is  particularly  crucial  for  tangential 
injection  that  carries  with  it  the  possibility  of  multiple  extrema,  overshoots, 
and  exponential  rather  than  algebraic  decay  of  the  wall  jet  velocity  profile  to 
the  freestream  levels. 

Referring  to  Fig.  44,  a  viscous  module  has  been  developed  to  described 
the  boundary  layers  in  Regions  (Tj)  ,  (?)  ,  (T),  (?),  and  the  blown  wall  jet  (?).  In 
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a  later  section,  prel  iminary  results  fora  module  for  region  (T)  will  be  dis¬ 
cussed.  Zones  (v)  and  (?)  are  in  the  vicinity  of  the  nose  and  are  especially 
important  to  treat  properly  in  the  small  disturbance  context  of  our  effort.  The 
correct  synthesis  and  treatment  of  the  small  disturbance  singularity  near  the 
stagnation  point  will  be  discussed  and  has  important  implications  regarding  the 
boundary  layer  development  in  (P)  and  (£^  .  Thus  far,  numerical  investigations 
suggest  that  the  precise  flow  details  in  (T)  and  (?)  are,  however,  not  decisive 
in  their  influence  on  the  global  structure  of  the  field,  in  particular,  the  role 
of  blowing  in  delaying  transonic  separation. 

On  this  basis,  a  rough  model  for  the  nose  region  combining  the  essential 
features  of  incompressible  stagnation  point  flow,  and  further  downstream,  small 
disturbance  singular  behavior,  was  implemented.  More  precise  simulations 
involving  asymptotic  and  computational  matching  of  the  aforementioned  singular 
solution  as  an  inner  limit  of  the  Karman  Guderley  outer  flow  and  a  full  potential 
flow  over  a  parabolic  nose  appropriate  to  finite  curvature  airfoils  have  been 
developed  by  K.  Kusonose  in  a  Ph.D.  dissertation  supervised  by  J.D.  Cole.^ 

This  solution  is  appropriate  for  a  more  accurate  description  of  the  inviscid  edge 
conditions  of  the  boundary  layer.  A  further  refinement  is  to  incorporate  these 
boundary  conditions  in  a  compressible  viscous  stagnation  point  flow  matching  with 
the  boundary  layer. 


4.2  Implementation  and  Results 

By  using  the  previously  enumerated  concepts,  a  laminar  boundary  layer  module 
has  been  constructed  based  on  a  generalization  to  tangential  blowing  of  the  box 
method  of  Ref.  34.  Results  have  been  computed  by  coupling  this  element  to  the 
inviscid  framework  previously  discussed  to  provide  a  blown  boundary  layer 
algorithm.  Formulational  details  are  given  in  the  appendices.  Specifically, 

Appendix  A  delineates  a  laminar  framework,  in  which  the  normal  coordinate  is 
written  as  a  Blasius  like  variable  which  Is  advantageous  in  resolving  the  steep 
gradients  and  treating  the  singular  flow  near  the  nose.  Appendix  B  describes 
special  procedures  Implemented  at  the  nose  that  are  required  to  properly  interface 
with  the  small  disturbance  inviscid  part  of  the  flow.  In  Appendix  C,  a  derivation 
of  jump  conditions  across  blown  wakes  according  to  matched  asymptotic  expansions 
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is  given.  These  relations  are  of  importance  to  future  efforts  on  characterizing 
viscous  modifications  of  the  external  flow. 

Briefly,  for  the  results  to  be  discussed,  a  finite  difference  algorithm 
based  on  the  compressible  boundary  layer  equations  was  used.  Furthermore,  a  unity 
Prandtl  number,  Chapman  gas  with  constant  density  viscosity  product  across  the 
boundary  layer,  and  Insulated  wall  was  assumed.  With  these  assumptions,  an 
integral  of  the  energy  equation  Implies  that  the  stagnation  enthalpy  in  units  of 
its  value  at  the  boundary  layer  edge,  E,  is  a  constant  =  1  across  this  layer.  The 
momentun  equation  is  thereby  decoupled.  Thermodynamic  variables  such  as  the 
density  in  the  layer  can  be  calculated  from  the  boundary  layer  assumption  of 
constant  pressure  across  the  layer  and  E  =  1.  These  observations  lead  to  the 
following  momentum  equation  for  the  normalized  velocity  f'  in  units  of  its 
boundary  layer  edge  value: 


(bf")'  +  nijff"  -  n^f'2  +  m^c  =  m1Q  |f' 


where  the  coefficients  (defined  in  Appendix  B)  are  related  to  the  external 
pressure  gradient,  primes  refer  to  partial  differentiation  with  respect  to  c ,  a 
scaled  Blaslus  coordinate  In  a  direction  perpendicular  to  the  parallels  to  the 
airfoil.  The  quantity  x  refers  to  the  coordinate  along  the  parallels  as  shown  in 
Fig.  44. 

Figure  45  depicts  a  NACA  0012  airfoil,  with  a  smoothed  small 
disturbance  SLOR  chordwise  pressure  distribution  on  its  upper  surface  for  =  0.7 
and  «  =  3°  as  a  basis  of  illustration;  typical  results,  shown  in  Fig.  46, 
indicate  the  streamwise  evolution  of  velocity  profiles  upstream  of  the  slot  and 
uncorrected  for  viscous  interaction  effects  on  the  external  flow  field.  Despite 
the  rather  severe  adverse  pressure  gradients,  particularly  those  associated  with 
the  nose  singularity,  and  their  significant  influence  on  the  source  term  In  the 
momentum  equations  as  the  coefficients  of  the  reduced  form  in  the  Blasius 
variables,  the  box  scheme  is  robust  enough  to  treat  such  variations.  For  the  case 
at  hand,  the  loss  In  fullness  In  the  profile  resulting  from  the  adverse  pressure 
gradient  is  evident.  Associated  decreasing  wall  shear  stress  is  also  apparent. 
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Fig.  45  Pressure  distribution  on  upper  surface  of  NACA  0012  airfoil,  M  =  0.7, 
a  =  3 
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Regarding  the  influence  of  slot  blowing.  Fig.  47  shows  evolution  for  the 
same  airfoil  of  the  profiles  downstream  of  a  slot  located  at  x*  =  0.2  in  which 
(x*,y*)  refer  hereinafter  to  Cartesian  coordinates  erected  at  the  mean  camber  line 
in  the  usual  way.  For  the  examples  selected,  an  initial  parabolic  slot  profile 
was  utilized,  in  which  the  velocity  function  at  the  slot  station  x*  =  x*  is  given 
by 


-  AU  -cdk  .  0  <  c*  <  4d 


» 


C  =  4  +  cd 


(2a) 


and  for  the  range  of  r*  above  the  slot  lip, 

)  =  f,(x*-,c*)  ,  cd  <  c*  <  •  (2b) 

where  the  solution  which  had  been  marched  from  upstream  to  the  slot  location  is 
assumed  to  have  a  continuous  velocity  component  at  the  slot  (excluding  the 
pathological  case  of  shocks  at  the  slot  location). 

From  (2a), 

3  2 

f  =  .  +  Acd  4-  •  0  <  4*  <  4d  (3a) 

A4d 

f(4d)  — 5“  •  (3b) 

For  more  general  slot  profiles  in  which 


f'(xs,c  )  =  F'(c  )  ,  0  <  q  <  5d  , 


the  blowing  coefficient  =  Cj  =  — ^ —  »  where  J  =  momentum  flux/unit  span, 
c  =  chord  is  given  by  P®LI  c 
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0  *:  * 
F'^U  )d? 


(4) 


* 

where  Re^  is  the  freestream  Reynolds  number  based  on  the  chord,  and  x$L  is  the 
streamwise  position  of  the  slot  in  units  of  the  chord.  Actually,  (4)  is  an 
approximation  to  within  terms  of  O(s^).  In  (4),  the  quantity  is  given  by: 


where  d  is  the  dimensional  slot  height  and  p  is  the  mean  density  across  the  slot. 


For  the  example  indicated  in  Fig.  45,  the  slot  was  located  at  20% 
chord,  very  close  to  the  onset  of  separation,  which  occurred  at  27.6%  chord  for 
this  case.  It  should  also  be  noted  that  the  location  of  the  shock  which  had 
been  slightly  smoothed  for  purposes  of  initial  checkout  of  the  algorithm  (which 
can  handle  non-smooth  cases)  is  at  about  23%  chord.  The  characteristic  diffu¬ 
sion  of  the  profile  associated  with  mixing  of  the  jet  is  evident  in  the  figure. 
More  significant  and  not  clearly  indicated  (but  shown  later),  is  the  fact  the  - 
the  separation  point  has  now  been  moved  downstream  to  79%  chord.  Other  calcula¬ 
tions  show  that  with  modest  further  increases  in  ^d,  separation  can  be  complete- 
ly  eliminated.  (Note  in  this  context  that  the  peak  velocity  for  (2a)  is  Acd/4.) 

In  Ref.  36  the  corresponding  profile  development  for  an  "underblown" 
case,  where  the  peak  is  less  than  the  freestream  value,  is  discussed.  Relaxation 
to  a  conventional  profile  occurs  as  previously.  However,  this  blowing 
configuration  actually  results  ir.  premature  separation,  compared  to  the  case 
of  no  blowing.  This  is  indicated  in  Fig.  48  where  the  effect  of  blowing  on  the 
separation  point  location  is  shown.  Respite  the  initially  higher  shear  stress 
at  the  slot  in  the  underblown  case  ^  =  1,  A  =  2,  the  higher  vorticity  diffusion 
and  lower  overall  momentum  in  the  layer  leads  to  earlier  separation  which  can  be 
seen  as  the  leftmost  solid  circle  in  Fig.  48.  This  level  should  be  related  to 
the  unblown  result  shown  in  dashed  lines  in  the  same  plot.  Moreover,  the  case 
Cd  =  2.65,  A  =  1,  moves  the  separation  point  substantially  downstream  to  almost 
80%  chord.  Other  cases  are  shown  in  Fig.  48,  such  as  the  one  corresponding  to 
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the  slot  position  at  approximately  11%  chord.  This  demonstrates  the  dramatic 
role  of  the  slot  location  and  upstream  boundary  layer  thickness  in  delaying 
separation;  where  for  the  same  slot  height,  the  upstream  slot  location  gives  an 
almost  trivial  downstream  movement  of  the  separation  location,  in  contrast  to 
the  potent  effect  of  the  downstream  slot  position. 

The  results  of  Fig.  48  can  be  replotted  as  in  Fig.  49  to  show  the 

trends  as  a  function  of  a  momentum  flux  parameter  which  is  proportional  to  Cj. 

By  Fq.  (3),  Cj  =  (x^/Re)1^2  A2r,J/30  for  a  parabolic  profile.  For  the  case  at 

hand,  according  to  Fig.  39,  an  A2^  =  200  represents  a  value  sufficient  to  drive 

the  separation  point  to  the  trailing  edge.  For  ReM  =  10^  and  p/p  =  1,  this 

"00 

implies  that  Cj  «  3  x  10”  ,  a  value  that  is  surprisingly  modest  in  view  of  the 
ease  in  which  laminar  boundary  layers  separate  and  the  typical  turbulent  values  of 
Cj  associated  with  transonic  experiments.  These  values  which  are  betweeen  0.01 
and  0.1  are  selected  to  be  large  enough  to  provide  significant  supercirculation 
enhancements.  Returning  to  the  example  calculated,  Eq.  (4)  indicates  that  the 
slot  height  d/c  =  0.9  x  10"3,  corresponding  to  a  «  2.5.  For  the  limited 
number  of  cases  run,  there  is  a  suggestion  in  Fig.  49  that  the  curves  of  Fig.  48 
collapse  to  a  single  universal  band  of  results.  In  view  of  the  roles  of  the 
pressure  gradient,  and  the  multiple  extrema  in  the  velocity  profile,  this 
assertion  must  be  regarded  as  tentative  at  best,  pending  more  extensive 
parametric  studies.  What  is  significant,  however,  is  that  for  the  first  time, 
the  delay  effect  of  tangential  blowing  on  natural  and  shock  induced  separation 
over  transonic  airfoils  has  been  inexpensively  quantified. 
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MOMENTUM  FLUX  PARAMETER  =  A2^ 


Fig.  49  Reduced  plot  of  separation  location  x*  as  a  function  of  momentum  flux 
parameter  A2 5, 5.  s 
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5.0  WAKE  STUDIES 


I 


I 


5.1  Formulation 

To  obtain  the  wake  structure,  Eq.  (1)  of  Section  4.2  is  solved  by  using 
"initial"  conditions  at  the  trailing  edge  from  the  boundary  layer  solutions 
obtained  with  our  parabolic  solver  upstream.  In  contrast  to  the  no  slip  and 
asymptotic  conditions  required  to  complete  the  formulation  over  the  airfoil,  the 
wake  requires  the  following  conditions  to  be  fulfilled: 


f  ■  (x,±oo)  = 
f (x,0)  =  0 


(la) 

(lb) 


where  x=l  is  the  location  of  the  trailing  edge. 


Note  that  the  zero  argument  in  (lb)  is  a  consequence  of  Prandtl's 
transposition  theorem,  where  the  boundary  condition  on  the  dividing  streamline 
can  be  replaced  by  those  at  any  other  convenient  location,  e.g.,  ^=0. 

With  the  Box  Scheme,  Eq.  (1)  gives  a  three  point  boundary  value  problem 
for  an  inhomogeneous  ordinary  differential  equation  formulation  based  on  the 
parabolic  partial  differential  Eq.  (1)  of  Section  4.2.  With  pivoting  strategies 
to  avoid  singular  matrices  such  as  those  developed  by  H.B.  Keller,  the  iterative 
inversion  of  the  associated  block  matrices  can  be  achieved  to  satisfy  (lb) 
directly.  For  coding  expendiency,  we  have,  however,  elected  to  replace  direct 
satisfaction  of  (lb)  by  an  iterative  approach,  in  which  we  introduce  an 
eigenvalue  y  Tor  which 


f(x,-»)  »  y(x) 


In  our  method,  the  quantity  y  is  determined  iteratively  so  that  (lb)  is 
sati sfied. 
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Numerical  Details 

On  the  upper  and  lower  surfaces  of  the  airfoil  we  solve  the  momentum 
equation  (Eq.  (1)  of  cection  4.2)  with  boundary  conditions 


fu(x.O)  = 


3ftt(x  ,0) 


u .  „  3VX'*> 


3n 


-  0, 


3n 


=  1 


and  f  (x,0)  =  sf  (x,0)/3n  =  0,  3f  (x,-»)/3n  =  1  on  the  upper  and  lower  surfaces, 

XX  X 


respectively.  For  x  =  x-jy,  we  form  the  vector 


where  ?u  and  t  are  the  discrete  solution  vectors  obtained  from  solving  the 
momentum  equation.  For  x  >  XyE,  we  continue  to  solve  this  equation  with  initial 
conditions  ^(x )  and  boundary  conditions  F(x,-»)  =  y(x)  as  well  as 

-{£«*. -i-i  . 


Here,  y(x)  is  selected  in  a  Newton-like  iterative  scheme  to  ensure  that 
F(x,0)  =  0,  i.e., 

Jn+D  -  v(»)  F(»,0, ■,(">] 

^(x.0.>,S 

where  the  superscript  is  the  iteration  counter.  In  practice,  we  approximate  the 
derivative  by 

3F(x  0  ,  F(x,0,y1">)  ■  FU,0 

3rlx,0,y  >  T(nj  :  Yt"-i) 

One  computational  difficulty  encountered  was  that  the  initial  condi¬ 
tions  were  not  smooth  enough  to  guarantee  convergence.  Accordingly,  a 

spline  least  square  smoothing  algorithm  was  employed  to  obtain  smoother  starting 
conditions. 
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5.2  Parametric  Studies 

Using  the  formulation  and  computational  procedure  described  in  the 
previous  section,  we  performed  a  series  of  studies  of  the  wake  structure 
behind  a  tangentially  blown  wake.  A  blown  NACA  0012  airfoil  at  MTO  =  0.7 
and  a  =  3°  was  investigated,  in  which  tangential  blowing  was  introduced  on  the 
lower  as  well  as  the  upper  surface  at  slots  located  at  21%  chord.  The  pressure 
distribution  on  the  latter  is  shown  in  Fig.  45.  Blowing  on  both  surfaces  was 
required  to  maintain  attached  flow  to  the  trailing  edge  in  view  of  the  adverse 
pressure  gradient  associated  with  the  (smoothed)  shock  recorapression  indicated 
in  Fig.  45. 

In  Fig.  50,  a  typical  streamwise  evolution  of  the  velocity  profile  u/ue 
downstream,  including  the  trailing  edge,  is  shown,  where  ug  =  the  external 
velocity  =>  U  the  freestream  value.  The  slot  height  is  expressed  in  terms  of 
the  Blasius  variable.  If  ue  is  the  velocity  at  the  boundary  layer  edge,  the 
peak  velocity  parameter  A  for  a  parabolic  velocity  profile  at  the  slot  where 

A  “  j4  ^max  at  slotj^d 

was  assigned  the  value  3.  The  other  parameters  for  this  case  are  indicated  in 
the  caption  of  this  figure.  Diffusion  to  the  symmetrical  bell  shaped  profile  is 
evident  in  the  figure.  Actually,  this  shape  was  achieved  as  early  as  x  =  5.96 
for  case  1  (not  shown).  In  spite  of  the  rapidity  of  this  relaxation  in  view  of 
the  considerable  initial  asymmetry,  the  viscous  interaction  process  could  be 
sensitive  to  the  details  of  streamwise  profile  evolution  in  the  transition 
phases. 

In  Fig.  51,  the  peak  velocity  (u/ue)max  as  well  as  displacement 
thickness  6  and  momentum  thickness  o,  in  units  of  c/Zffe^  where  c  =  airfoil  chord 
and  Rem  =  freestream  Reynolds  number,  are  shown.  By  contrast  to  the  classical 
submerged  jet  studied  by  Schl ichting^  and  Bickley^  for  incompressible  flows, 
the  momentum  thickness  in  this  coflowing  case  is  strongly  affected  by  the 
external  pressure  gradient  and  has  a  significant  decay  for  the  first  3  chord 
lengths  downstream.  Thereafter,  this  parameter  remains  almost  constant  as  in 
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Fig.  50  Streamwise  evolution  from  trailing  edge  of  wake  velocity  profiles  for 
NACA  0012  airfoil,  =  0.7,  o  *  3°,  blown  tangentially  on  upper  and 
lower  surfaces  at  slot  position  Xc^  of  21%  chord,  peak  velocity 
parameter  A  =  3,  normalized  slot  height  ^  =  2.5. 
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the  submerged  case,  in  spite  of  the  external  stream.  (As  an  approximation  to 
obtain  these  results,  the  pressure  coefficient  Cp  at  the  wake  edge  calculated 
with  an  SLOR  solution*  has  been  extrapolated  downstream  by  using  a  monomial 
fit  Cp  =  Ax^.)  Correspondingly,  the  displacement  thickness  relaxes  to  a  similar 
invariance  more  quickly,  as  indicated  in  Fig.  51. 

It  was  of  interest  to  compare  the  bell  shape  profile  to  an  asymptotic 
flat  plate  wake  solution.  Strictly  speaking,  this  comparison  is  only  of  valid¬ 
ity  when  the  reduced  perturbation  velocity  is  small.  Although  this  is  not  the 
case  at  the  streamwise  stations  under  consideration,  we  study  here  the  extended 
validity  of  this  behavior,  from  the  viewpoint  of  experimental  and  theoretical 
interest. 

Denoting  the  perturbation  velocity  (u  -  IJ ) / U  as  u',  where  U  is  the 
freestream,  we  define  g'  and  e  such  that 


f'  =  1  +  u'  =  n  eg  ,  (e  0) 
f  =  n  '  eg 


The  momentum  Eq.  (1)  of  bection  4.2  specialized  to  this  laminar  case  which  is 


* 1 r -  ■  *  (f ' ■  w 


t"  |Ii 

3X ' 


(2) 


where  the  coefficients  in  the  original  equation  have  taken  on  limiting  values  in 
the  asymptotic  nearly  zero  pressure  gradient  flow  reduces  to  the  linearized  form 


+  7  9"  ’  x  |x"  =  0 


f '  =  i-1 


(3a) 


with  the  boundary  conditions 


9(0)  =  g'(«)  =  g"(0)  =  0 


(3b) 


♦Assistance  of  V.  Shankar  in  obtaining  these  results  is  gratefully  acknowledged. 
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On  the  basis  of  the  invariance  of  e  and,  correspondingly,  the  momentum  flux  J, 
for  large  x  where 

"  2 

J  =  /  pu  dy  =  const.  ,  (4) 

-oo 

with  y  a  Cartesian  coordinate  normal  to  the  freestream,  p  =  density,  the 
quantity  g'  has  the  following  similitude 

9*  =  G(n)  •  (5) 


This  leads  to  the  well  known  solution  of  (3)  given  by 

a.  Ce-2/< 


(6) 


The  effect  of  the  finite  size  of  the  perturbation  velocity  u'  is  shown  in 

Fig.  52,  where  fhe  scaling  of  the  peak  velocity  u't,  =  eg  (x,0)  is  indicated. 

Pe°x  n  q  3Qin 

By  contrast  to  (6),  the  numerical  results  suggest  that  u'eak  -  x  *  .  If, 

however,  C  in  (6)  is  determined  by  matching  with  the  numerical  solution  such  as 
shown  in  Fig.  53,  or  alternatively,  u'  is  normalized  to  the  peak  velocity  at 
each  x  station,  a  strikingly  good  agreement  between  Eq.  (6)  and  the  computa¬ 
tionally  obtained  wake  profile  is  obtained,  as  indicated  for  the  conditions 
delineated  in  the  figure  for  the  blown  NACA  0CT2  airfoil. 

Regarding  these  facts,  we  have  also  determined  that  Upg-^  •-  x-8**885  as 
x  *  «  for  the  case  at  hand,  in  contrast  to  the  Bickley  submerged  jet  solution 
where  Upe,t.  x"^8.  The  slower  rate  of  decay  of  the  present  coflowing  jet 
seems  consistent  with  the  reduced  action  of  viscous  dissipation  associated  with 
the  lower  relative  velocity  between  the  jet  and  external  flow.  It  should  also 

OO 

be  noted  that  this  trend  differs  from  that  discussed  in  Birkhoff's  book 
regarding  wakes  with  "hydrodynamic  self  propulsion"  associated  with  self  pro¬ 
pelled  objects  such  as  boats  or  airplanes.  For  these  cases,  if  no  acceleration 
occurs,  the  thrust,  T,  is  equal  to  the  drag,  D,  leading  to 

00 

j  =  J  u'dy  =  T  -  D  -  0  ,  (7) 
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Fig.  52  Comparison  of  peak  wake  perturbation  velocity  (from  freestream) 

downstream  decay  against  flat  plate  similarity  solution  for  case  of 
Figures  50  and  51 . 


lengths  downstream 


5*Va» 


i 

i  ■*- 


I 


I 


SC5055.21FR 


For  such  cases,  it  can  be  shown  that 
00  ? 

M  =  /  y  u'dy  ,  invariant  with  x,  (8) 

-oo 

The  invariance  (8)  leads  to  a  different  similarity  than  for  the  Bickley 
jet.  For  self  propelled  wakes,  Upg^  -  The  weaker  decay  experienced  in 

the  present  studies  is  associated  with  the  lack  of  similarity  arising  from  the 
presence  of  two  characteristic  velocities,  i.e.,  that  from  the  jet  momentum,  and 
the  other  corresponding  to  the  freestreain.  The  latter  boundary  condition  and 
(4)  Impose  conflicting  requirements  negating  similitude.  The  other  reason  for 
deviation  from  the  self  propelled  case  is  the  fact  that  T  >  D  in  the  flow  being 
computed  here. 

In  future  studies,  two  regimes  will  be  considered  appropriate  to  the 
blown  wakes  considered.  For  T  »  D,  the  blown  wake  can  be  considered  as  a 
linear  nonsimilar  perturbation  on  the  Bickley  submerged  jet  similar  solution  for 
sufficiently  small  x.  For  x  ■+  »,  the  wake  behaves  like  a  classical  flat  plate 
structure  in  which  T  -  D  is  conserved. 

Besides  their  theoretical  significance,  these  considerations  are  of 
importance  in  assessing  the  net  thrust  of  propulsive  wings  in  experimental 
simulations.  In  addition,  the  asymptotic  wake  characteristics  discussed  here 
will  be  employed  in  modeling  the  far  field  for  the  blown  viscous  interaction 
problem  in  subsequent  studies. 


i 
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6.0  CONCLUSIONS 


During  the  course  of  the  research  program,  a  number  of  areas  received 
attention  in  connection  with  the  study  of  tangentially  blown  airfoils.  Par¬ 
ticular  emphasis  has  been  placed  on  understanding  the  structure  of  compressible 
wall  jets  occurring  in  such  flow  fields.  Modern  relaxation  methods  have  pro¬ 
vided  an  efficient  means  of  analyzing  the  diversity  of  wave  patterns  arising  in 
these  jets.  For  submerged  transonic  jets,  with  subsonic  conditions  infinitely 
far  downstream  and  a  specified  exit  velocity  profile,  the  following  specific 
findings  arise  associated  with  unchoked  flow  over  a  series  of  wall  shapes. 


'  I 

I 


•  The  Kutta  condition  on  the  nozzle  rim  can  be  satisfied  merely  by 
requiring  streamwise  continuity  of  the  potential  across  the  rim  as 
a  trailing  edge  for  the  flow.  With  the  use  of  this  procedure,  the 
numerical  algorithm  correctly  tracks  a  local  singular  solution 
which  has  a  square  root  zero  and  is  locally  harmonic  in  scaled 
variables,  providing  the  rim  is  not  transonic. 

•  As  compared  to  unconfined  flows,  the  slip  line  boundaries  create  a 
rapid  decay  of  the  disturbances.  The  functional  form  of  the  far 
field  perturbation  potential  is  an  exponential  damped  sine,  similar 
to  that  encountered  for  incompressible  flow  but  different  in  that 
its  amplitude  interacts  nonlinearly  with  the  near  field. 

•  Analysis  of  the  free  jet  case  indicates  subcritical  monotone 
streamwise  variations  of  the  pressure  for  a  subsonic  jet  exit,  as 
in  linear  subsonic  flow. 

•  Acceleration  of  the  wall  jet  to  criticality  over  convex  walls  is 
accomplished  by  stream  tube  contractions  and  throats  induced  by 
upstream  influence  of  the  turning. 
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•  The  upper  slip  line  of  the  wall  jet  becomes  asymptotically 
parallel  to  that  emanating  from  the  trailing  edge  infinitely  far 
downstream. 

For  this  phase  of  the  effort,  as  well  as  that  dealing  with  coflowing  jets,  the 
inverse  methods  reported  in  Refs.  27-29  appear  to  be  an  attractive  means  to 
reduce  wave  drag  and  enhance  thrust  recovery  by  elimination  of  shock  waves  by 
suitable  wall  contouring. 

With  regard  to  the  investigation  of  coflowing  wall  jets,  selection 
rules  have  been  identified  for  the  existence  of  throats  in  the  jets  on  the  basis 
of  an  integral  theorem  based  on  small  disturbance  formulation  of  a  unit  problem 
involving  a  flat  plate  cowl.  The  latter  has  important  features  that  are  of 
relevance  to  more  complicated  cowl  geometries.  For  a  slightly  subsonic  external 
flow,  with  an  embedded  jet  at  nearly  the  same  velocity  and  settling  pressure, 
approximately  sixteen  cases  have  been  enumerated  depending  on  the  range  of  jet 
similarity  parameter.  Computational  results  have  validated  some  of  these  cases. 
The  information  provided  in  this  phase  illuminates  factors  influencing  the  pene¬ 
tration  length  of  initially  supersonic  jets  and  the  scales  affecting  overexpan¬ 
sion  spikes  evident  ,-n  experiments  over  blown  surfaces  of  propulsive  wings. 

In  the  portion  of  the  study  dealing  with  the  complete  problem  of  a  tan¬ 
gentially  blown  airfoil,  asymptotic  and  computational  simulations  have  received 
attention  in  incompressible  and  transonic  regimes.  These  analyses  have  utilized 
models  of  the  jet  that  are  accurate  in  regions  away  from  the  nozzle  exits  and  that 
are  capable  of  later  integration  of  the  results  of  the  aforementioned  basic 
studies  of  coflowing  wall  jets  appropriate  to  a  local  description  of  the  exit 
flows.  For  the  more  globally  oriented  inviscid  blown  airfoil  models,  the 
analytical  and  computational  models  Indicate  that: 

•  In  the  thin  jet  small  deflection  approximation,  the  pressure  jumps 
associated  with  the  Spence  theory  prevail  even  if  the  flow  is 
rotational  and  compressible. 
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•  The  asymptotic  developments  provided  allow  further  systematic 
refinements. 

•  Effects  associated  with  initial  skewness  and  vorticity 
inaccessible  to  other  theories  can  be  assessed. 

•  For  finite  trailing  edge  angles,  the  dividing  streamline  leaves 
tangent  to  the  higher  stagnation  pressure  side  if  upper  and  lower 
stagnation  pressures  are  not  equal  and  at  the  trailing  edge 
bisector  if  they  are. 

•  Computational  results  obtained  for  transonic  USB  configurations 
indicate  significant  enhancements  of  lifting  pressures  associated 
with  blowing. 

•  Comparisons  with  experiment  indicate  the  need  for  refinements 
incorporating  wave  interaction  phenomena  near  the  jet  exit  as  well 
as  viscous  interaction  processes  in  the  downstream  portion  of  the 
wall  jet. 

Initial  results  characterizing  viscous  effects  using  finite  difference  laminar 
descriptions  of  the  flow  fields  demonstrate  clearly  that: 

•  Substantial  downstream  movements  in  the  shock-induced  separation 
point  are  achievable  with  application  of  tangential  blowing.  By 
computational  schemes,  these  delays  can  be  inexpensively 
quantified  as  compared  to  experimental  methods. 

•  The  propulsive  wake  at  great  distances  downstream  gives  a  nearly 
Gaussian  normalized  velocity  profile  in  spite  of  a  lack  of  self 
similarity. which  occurs  in  purely  submerged  or  self  propelled 
cases.  This  result  is  of  significance  for  assessment  of  net' thrust 
in  experimental  simulations  and  theoretical  far-field  modeling. 
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With  the  viscous  module  and  the  appropriate  Iterative  coupling  algor¬ 
ithm  to  the  external  flow  to  be  implemented  in  the  near  future,  optimization 
between  separation  suppression,  wave  drag  minimization,  and  supercirculation 
control  will  be  possible.  It  is  envisioned  that  the  design  techniques  contained 
in  Refs.  27-29  will  augment  this  capability  by  providing  methods  to  modulate 
shock  formation  In  concert  with  the  blowing  effects. 

On  a  more  ambitious  scale,  which  will  depend  on  further  developments  in 
current  parabolic  marching  schemes  and  improvement  and  application  of  multiple 
deck  asymptotic  models,  future  effort  should  consider  blowing  effects  on: 

•  Strong  shock  boundary  layer  interactions  with  limited  separation 
near  the  foot  of  the  Incident  normal  shock  and  trailing  edge 

•  Blowing  effect  on  near  wake  curvature  modification  of  viscous  lift 
at  trailing  edge  and  cowl  lip 

•  Impact  of  shock  waves  in  the  jet  on  wall  jet  transition  and 
separation. 
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APPENDIX  A 

DETAILS  OF  VISCOUS  FORMULATION 


Governing  Equations 

In  what  follows,  the  approach  to  treat  the  viscous  regions  of  the  flow 
field  will  be  described.  In  major  respects,  the  formulational  aspects  given 
here  resemble  those  given  in  Ref.  3$.  However,  certain  important  details  asso¬ 
ciated  with  the  blowing  application  require  special  treatment.  For  a  self  con¬ 
tained  account,  therefore,  some  of  the  unblown  formulation  is  repeated  here. 
Referring  to  the  curvilinear  axis  system  shown  in  Fig.  Al,  the  boundary  layer 
equations  (which  include  laminar  and  turbulent  cases)  to  be  utilized  to 
calculate  viscous  zones  are  as  follows: 


Continuity 

lx  <»“>  *  !y  C^l>  *  0 


(la) 


x-Momentum 


where  as  a  typical  average,  ^v  =  pv  +  p1 v’1 ,  and  as  usual,  the  primes  denote 
fluctuations  about  temporal  averages,*  Pr  is  the  Prandtl  number,  p  =  density, 
u,v  are  velocity  components  in  the  x  and  y  direction,  respectively,  H  =  total 


*The  appropriate  eddy  viscosity  representations  for  the  Reynolds  stresses  have 
not  been  considered  in  this  phase  of  the  effort,  but  will  receive  attention 
in  later  aspects  concerning  marching  of  the  boundary  layer  downstream  of 
impingement  points  of  shocks  from  the  external  flow. 
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JET  AND  WAKE 


Fig.  A1  Upper  surface  blown  airfoil  showing  viscous  zones. 


146 


Rockwell  International 

Sctenc*  Center 


SC5055.21FR 

enthalpy,  \i  =  dynamic  viscosity,  and  h^  =  metric.  If  in  the  x*,y*  Cartesian 
coordinate  system  shown  in  Fig.  Al,  the  surface  of  the  airfoil  Is  given  as 

y  =  F(x*)  , 

then  the  metric  h-  is  given  by 

hj  =  /r+rrwrp  . 

Boundary  Conditions 

The  appropriate  conditions  for  the  solution  of  (1)  consist  of  surface 
constraints,  initial  conditions,  and  matching  considerations  to  the  outer  flow. 
The  surface  conditions  are  given  as  follows: 

Surface  (No  Slip  and  Adiabatic  Wall) 


u(x,0)  =  v(x,0)  =  0  (2a) 

~  (x  ,0)  =  0  .  (2b) 

Initial  Conditions 

In  Appendix  B,  a  local  nose  solution  is  considered  which  Is  utilized  to 
start  the  solution  In  a  manner  appropriate  to  the  parabolic  initial  boundary 
value  problem  for  Eqs.  (1). 

At  the  jet  exit,  x  =  x^,  referring  to  Fig.  A2;  an  initial  profile  is 
specified  on  the  blown  section  AB.  The  data  on  BC  is  "inherited"  from  history 
of  the  flow  field  which  has  been  marched  to  this  point  from  the  boundary  layer 
and  external  flow  through  iterations  between  successive  sweeps  of  some  sort  of 
successive  line  relaxation  method  or  other  inviscid  solver.  A  key  point  is  the 
following:  The  solution  of  (1),  subject  to  the  initial  conditions  in  ABC,  char¬ 
acterizes  the  mixing  of  the  flow  exhausting  from  the  nozzle  and  the  boundary 
layer  along  the  upper  surface  of  the  airfoil  (ORC  in  Fig.  Al).  A  similar  treat- 
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ment  has  been  used  to  handle  the  mixing  of  the  upper  and  lower  streams  in  the 
wake  region  downstream  of  the  trailing  edge  and  to  establish  the  jump  conditions 
for  the  pressure.  These  will  be  a  viscous  modification  of  the  Spence  jump  con¬ 
ditions  developed  from  concepts  of  Preston  and  implemented  by  us  (Refs.  5-7). 

Matching  with  Outer  Flow 

The  appropriate  boundary  conditions  at  the  edge  of  the  viscous  region 

are: 


lim  u(x,y) 
y  +  °° 

-  ue<x> 

(3a) 

lim  v(x,y) 
y  + 

»  ve(x) 

(3b) 

lim  p(x,y) 
y  +  * 

=  pe(x)  , 

(3c) 

where  the  quantities  on  the  right  side  of  Eqs.  (3)  are  effectively  the  inner 
limit  of  the  outer  flow,  i.e.,  if  W(x,y)  is  one  typical  such  quantity,  then. 

W  (x)  »  lim  W(x,y).  (e  -  external.) 
y+« 

Transformed  Problem 

In  this  section,  we  generalize  the  procedure  given  in  Ref.  2  to  the 
blown  boundary  layer  problem  indicated  previously.  On  Introduction  of 
transformed  independent  variables 

/  ue  \  1/2  x 

x  =  x,  dc  =  ----- - )  pdy  ,  s.  =  /  h.dx  (4a) 

\pewesl/  1  0  1 

and  the  dependent  variables 

<4b> 

with 

2  /2 

*  =  (pe^eVl*  ’  (4c) 
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the  system  (1)  under  the  transformation  (4)  becomes 


(bf")'  +  n^ff"  +  m2f'2  -  mnc  =  m1Q  [f '  -  f" 


i  df  ^ii  3f 
3X  "  ax 


<“1E')'  +“2E'  +  “3'"lO  (f'  H-  E'  w) 
where  primes  denote  partial  differentiation  with  respect  to  5,  and 


peue 


(5a) 

(5b) 


=  eddy  viscosity  (=  0  for  laminar  flow) 


e  signifies  conditions  on  the  edge  of  the  boundary  layer 

^  =  ?  (E  ♦  iqu;  sr)  *  h  (/v?  ■  ?  (1  *  "z> +  I?  (/v? 

X 

s 1  -  J  h.dx 
0  1 

m  S1  3ue 
m2  hjU0  3x 


*10  =  h 


S1  3ue 

11  *  “e15!  \  "  1  t 

M1  =  PF  y  +  4  Pr^J  ’  n2  =  mlf  »  y3  =  c  ^  "  Pr  ^  f'f" 

Here,  dagger  superscripts  refer  to  turbulent  quantities.  Herein,  it  will  be 
assumed  for  convenience  that  the  Prandtl  number  is  unity. 
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Boundary  Conditions 

The  transformed  versions  of  (2)  and  (3)  are 
c  =  0  ,  f  =  f'  =  E'  =  0  (6a) 

U  »  O  »  1)  ,  f'  -  E'  -  1  •  (6b) 

In  accord  with  the  remarks  given  in  Section  4.1,  we  have  assumed  in  our  studies 
that  Pr  =  1,  C  =  1  (Chapman  gas).  An  integral  of  (5b)  satisfying  these 
conditions  is  E  =  1.  In  addition,  we  assume  throughout  that  ue  =  where 
subscript  »  denotes  freestream  conditions. 


Numerical  Procedures 

The  discretization  of  the  system  (5)  appropriate  to  treatment  of  the 
viscous  regions  will  be  Indicated.  These  equations  are  written  as  a  first  order 
system  by  introducing  the  new  independent  variables  u(x,c),  v(x,i;),  and  G(x,c). 
Thus,  Eqs.  (5)  can  be  written  as 


f*  »  U 

(7a) 

u'  =  V 

(7b) 

(bv)'  +  n^fv 

p 

+  (^(c-U  ) 

■  "lot11 

Hi  .  v 
ax 

Hi 

ax' 

(7c) 

E'  =  G 

(7d) 

(u^G)'  + 

+  ^  =  m10 

(u  —  - 

^  ax 

ell) 

ax-1 

• 

(7e) 

Using  Keller's  box  scheme  (7a),  (7b),  and  (7d)  become 
fn  -  fn 

j 

—  V/2 

n  J„n 

UJ  "  jzl  =  v" 
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(8c) 


where  we  have  averaged  about  the  point  (xn»Cj_i/2)*  Using  the  box  centered  at 
(xn_i /2 »Cj_i /2) »  Eq.  (7c)  is  discretized  as: 


+  (mi)n-W2fj.i/2vj.i/2  +  (m2^n-1/2^c"u2^j-1/2 


feWf  fiW 


II 


Rockwell  International 

Sci«nc«  C*nt«r 


m*:!  -  .  (^j‘:r  •  (*3>i" 

= - - + - f r - 

J  J 

,  \  '>,.n-i/2rn-l  x  rn“l  /*n  *n-l  1  L 

=  ^uj_l/2hj.l/2  +  bj.l/2  lTj_l  12  ~  T  j  _  1  /  2  J  jK  n 
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r  sn-1/2 


The  boundary  conditions  for  system  (8)  are 

fS  -  uS  - 0  u3  - 1  (9a> 

Gj  =  0  Ej  =  1  .  (9b) 

Note  that  (8c)  and  (8e)  form  a  linear  system  which  is  decoupled  from  (8abc). 
Therefore,  we  first  solve  the  nonlinear  system  (8abc)  with  boundary  conditions 
(9a).  Then,  we  solve  the  linear  system  (8de)  with  boundary  conditions  (9b). 
Initial  conditions’  are  obtained  from  considerations  in  Appendix  B  and  the  jet 
exit  conditions. 
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APPENDIX  B 

NOSE  SOLUTION  AND  SPECIAL  PROCEDURES 

The  approximate  uniformly  valid  solution  utilized  applies  to  a  wide 
class  of  airfoils  of  finite  radius  of  curvature.  Denoting  x*,y*  Cartesian 
coordinates  as  shown  in  Fig.  44,  this  family  of  airfoils  has  a  local 
representation  near  the  nose  given  by 

y*  =  fi^Zx*  +  a^x*  +  0(x*2)]  (1) 

where  aQ  and  a^  are  constants  and  6  is  the  thickness  ratio.  For  example,  the 
NACA  0012  has 


6  =  0.12  ,  aQ  =  1.4845  ,  a^  =  0.63  -  a 

where  a  Is  the  incidence,  and  all  lengths  are  units  of  the  chord. 

Based  on  the  previous  remarks,  a  uniformly  valid  representation  for  the 
edge  velocity  component  parallel  to  the  surface  ue  was  assumed  to  be 


where  the  constant  Cq  is  determined  by  patching  to  the  SLOR  solution  and  U  is 
the  freestream  velocity.  Note  that  as  x*  +  » 


V 


which  is  in  agreement  with  the  inner  expansion  of  the  Karman  Guderley  flow.  On 
the  other  hand,  as  x*  +  0 


FRECSXflNtt  PAG#  BlANK-NCT  FILMED 
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e  2 

-g-  *  Ax  +  Bx  +  •  ••  (Incompressible  stagnation  point  flow) 
x  =  arc  length 


with 


We  note  that  (1)  can  be  reverted  to  give 
*  i  M2 

ft*  =  IT  +  7?  4-r  +  (x3)  • 

6a0  6  a0 

From  (2)-(4),  the  pressure  gradient  parameter 
x  8up 

G<*>  -=  re»r 

appearing  in  the  rt^'s  defined  in  Appendix  A  can  be  evaluated  as 


(4) 


B 


G(x)  =  1  +  j  x  +  <?(x2) 


Furthermore,  the  quantity  c  appearing  in  Eq.  (5a)  is  obtained  as 
1  -  mf'2 


c  = 


T  -"m 


(5a) 


with 


m  = 


dr)2 


i  + 


(5b) 


In  our  procedure,  the  foregoing  relations  in  this  section  are  utilized  to  obtain 
a  boundary  layer  solution  from  the  stagnation  point,  xs,  to  a  downstream  station 
Xj.  For  x  >  Xj,  ue  Is  determined  from  the  SLOR  solution  and  the  boundary  layer 
equations  are  matched  using  the  data  at  Xj  as  initial  conditions. 
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APPENDIX  C 

JUMP  CONDITIONS  ACROSS  BOUNDARY  LAYERS  AND  WAKES 

In  this  section,  the  appropriate  jump  conditions  across  the  viscous 
regions  will  be  studied  by  using  second  order  boundary  layer  theory.  For  this 
purpose,  matched  asymptotic  expansion  procedures  will  be  employed  by  using 
intermediate  variables  to  effect  the  matching.  The  derivation  indicated  will  be 
for  incompressible  conditions.  However,  similar  results  are  envisaged  for  the 
compressible  case. 

In  Fig.  Cl,  a  coordinate  system  consisting  of  parallels  to  the  surface 
and  their  orthogonal  trajectories  is  indicated.  The  Navier-Stokes  equations  for 
incompressible  flow  are: 

Continuity:  v  •  q  =  0 

2 

Momentum:  -  §xu=-q*V$=--^-v  curl  u 

where  w  =  vorticity  =  curl  5  =  velocity  vector,  p  =  pressure,  and  \>  = 
kinematic  viscosity.  In  the  coordinates  indicated,  denoting  s  as  the  arc  length 
along  the  parallels  t  =  constant, 

ds1  =  differential  arc  length  along  surface  =  R(s)de 
R  =  radius  of  curvature,  o  =  angle  between  normals,  s  =  constant 
ds  =  differential  arc  length  along  arc  above  surface  =  (R+t)do. 

For  this  system,  an  elementary  displacement  vector  d#.  can  be  expressed 
in  terms  of  the  metrics  h^  and  as 

d^  =  h2ds?  +  h|dt2 


where 


hl  =  TT 


*(s)' 


h2  =  1 
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If  qt  and  qs  represent  the  velocity  components  in  the  t  and  s 
directions,  respectively,  then  the  vorticity  to  can  be  expressed  as 


where  £  is  normal  to  the  x,y  plane,  with  x  and  y  being  Cartesian  coordinates  as 
indicated  in  Fig.  Cl.  The  quantity  tu  is  then  given  by 


1  [3qt  3_  n  t  . 
“  ”  j  +  t  as  at  (  +  R(sT  qs) 
IT 


Moreover, 

_  - ,  _-j  -►  3u  f-  1  3u>  #• 

curl  „  =  __  r$  +  _ _  rt 
1  +  I 

where  and  are  unit  vectors  in  the  s  and  t  directions,  respectively. 

By  using  the  appropriate  relations  for  the  divergence  and  gradient  in  this 
curvilinear  coordinate  system  and 

3  x  <o  =  -  +  qt“*s  » 

we  obtain  the  following  component  form  of  the  Navier-Stokes  equations,  where  the 
kinematic  viscosity  has  been  normalized  to  be  proportional  to  the  reciprocal  of 
the  Reynolds  number,  e: 


Continuity 


3_ 

3S  +  3t 


l^1  +  ^qtl 
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Momentum 


,  ^ 
Hs  3S 

•  + 

(> *  *K 

3qs  + 

3 r + 

Vs 
TT“  : 

=  - 

3P 

3S 

^s 

3qt 

3Qt 

.  +  n 

<s 

3P 

1  +  i 

3S 

qt  at 

R  + 

T  =  " 

at 

1 

1  3<*t 

aqs 

q, 

s 

UJ  — 

>43S 

aTT 

'  ITT 

T 

• 

d\>  ( i  ,  ts  3o) 

as  '  <l  +  T?)  at 


_e _ _3w 

t  as 


1  + 


(lb) 


(lc) 


The  boundary  layer  equations  for  this  curvilinear  network  are  obtained 
using  the  following  limit  process  asymptotic  expressions  for  the  flow  variables: 


qs  =  u^s.t*)  +  fl  u2(s,t*)  +  ... 

\  -  'ft  v^s.t*)  +  eY2(s.t*)  + 

P  =  P^S)  +  /I  p2(s,t*)  + 
e  0  ,  s,t*  fixed 

Substitution  of  (2)  into  (1)  gives  the  following  hierarchies  on 
retention  of  the  indicated  orders  of  terms: 


(2a) 

(2b) 

(2c) 

(2d) 


3u,  au,  3P,(s)  a2u, 

^(1):  U1  ‘aiT  +  V1  at*  =  as~  +  at*? 


(3a) 


3PX 

3^ 


0  =  -  TTT  =  >  Pi 


PX(s) 


(3b) 


(/(ft) 


3U,  3U,  3U,  3U,  .*  au,  U,Vi 

II  .  .  4.  y  ..  .  ■  4-  y  - Z.  +  y  _ Z.  4.  ...  y  .  A  4.  ...  A. 

ul  as  u2  as  vi  at*  +  v2  at*  R  v2  at*  r 


ap2 

ls~ 


(4a) 
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u?  3Po 

r/d): 

(4b) 

3U.  3V, 

(5a) 

au9  av9  .  t*v. 

<?(^) :  TT  +  af7  +  af7  HH  =  0  * 

(5b) 

For  the  outer  inviscid  flow,  the  appropriate  asymptotic  expansions  are 

written  as 

qs  =  Ul(s,t)  +  /e  U2(s,t)  + 

(6a) 

qt  =  Vj(s,t)  +  /e  V2(s,t)  +  ••• 

(6b) 

P  =  Pj(s,t)  +  /e  P2(s»t)  +  ••• 

(6c) 

for 

e  +  0  ,  t,s  fixed 

(6d) 

On  substitution  of  (6)  into  (1)  and  retention  of  the 
following  equations  for  the  approximate  quantities  in 

indicated  orders,  the 
(6)  are  obtained: 

au,  t  3Ut  3P* 

u,  ^  (1  -aT 

(7a) 

“i  ,  aVi  ui  3,,i 

,  t  as  +  vl  at  R  3t 

1  4 1 

(7b) 

au,  t  3V,  v, 

jr +  +  ^  JT  "  TT  =  0 

(7c) 
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3Up  3U,  „  /  3U,  3V. 


V1U2  +  V2V1 


3P, 

W 


3V?  3V, 

"lir^iir.  v  !!i+¥  !!i.2Vi 
,  .  t  wi  st  *2  3t  n? 

A  R 


aP2 

W 


3U 

3S 


p  i  sVp  Vp 

- ♦  (' +  ir)  ir 4  ir  ■ 0  • 


(8a) 


(8b) 


(8c) 


The  appropriate  matching  between  the  inviscid  and  viscous  regions  is 
obtained  by  writing  the  expression  for  both  regions  in  intermediate  variables 
and  equating  li>a  terms.  The  appropriate  intermediate  variable  is  defined  as 

t  =  -4-y  .  where  e  «  n(e)  «  1  .  (9) 
n  ole) 


Thus : 


t* 


_  t 


t  =  nt  -*■  0 
n 


t*  =  t  +  a>  . 
ve  n 


Matching 

By  using  the  process  indicated,  the  inner  (viscous)  and  outer  inviscid 
representations  are  matching  schematically,  indicated  by  <*■  (asymptotic  equality) 
as  follows: 

P1(s,ntn)  +  /e  P2(s,ntn)  <--►  Pj(s)  +  /7  p2  s,^  t  +  •  •• 


or 
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3P  3P 

Pj(s,0)  +  tyt^  (s,0)  +  •••  +  /e  P2(s,0)  +  Ti  +  ••• 

p^(s)  +  /e  Tim  p2(s,t)  +  ••• 
t-*» 


which  for  matching  implies 


Px(s)  =  P^s.O) 


3P 


P2(s,t)  =  TT  ^s‘0^  +  P2^s*0^ 


(10a) 

(10b) 


l 

.£ 


and  for  t  -►  <*> 

U?(s,0) 

P2(s,t)  f"  fc  +  p2^s,°^  *  (10c) 

For  the  vertical  velocity  component,  we  have  analogously 

V^s.nt^)  +  /e  V2(s,ntri)  +  ...++  /I  V^s,-^  t^)  +  ev^s,-^  +  ••• 

3  V 

Vj(s,0)  +  nt)(  (s,0)  +  /e  V2(s,0)  +  *••■«-*■  /e  vi(s»7£  ^r|)  +  .  (H) 

For  compatibility  with  the  boundary  conditions,  we  assert  that  V^(s,0)  =  0. 
From  (11),  we  have  finally. 


(  9V,  \ 

(s,0)  =  lim  <  Vj(s,t)  =  t*  (s,0)  ( 

t*-x»  (  ' 

=  lim  jv^s.t)  +  t*  ■—  (s,0)  | 


ds* 


V2(s»0)  =  lim  -rj  /  (U,(s,0)  -  u1(s,T))dT  = 

t*-W°  0 


(12) 


6* 


displacement  thickness  =  ^  [Ug  -  u^(s,i)]dT  ,  Ug  =  Uj(s,0) 
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where  by  virtue  of  (5a)  the  following  relation  has  been  used  to  obtain  the 
displacement  thickness  integral  in  (12) 

t*  au, 

V1  =  /  IT  (s‘T>dT  * 


Pressure  Jump  Conditions  on  Wake  and  Boundary  Layer 

Referring  to  Fig.  C2,  we  consider  the  section  AB  of  the  wake.  The 
viscous  interaction  with  the  inviscid  portion  of  the  solution  will  be  effected 
by  pressure  jump  conditions  related  to  the  momentum  defect  in  the  viscous 
layers.  This  defect  plays  the  role  of  the  classical  Spence  blowing  coefficient 
utilized  in  our  jet  flap  algorithm.  To  see  this,  we  recall  from  previously 
derived  relations  that 

3P 

P2(s,t*)  +  (s,0)t*  +  P2(s,tt-)  ,  t*  +  *> 

3p 

P2(s*t*)  (s.Ojt*  +  P2(s,0-)  ,  t*  +  -co 

3Po  U?(s,t*) 

aS*  <*•**>  =~W~  ‘ 

On  integration: 

t* 

P2(s,t*)  -  p2(s,0)  =  j  /  u|(s,t*)dt*  . 

A1  so, 

3P 

-jjA  (s»0+)  =  "ft  U|(s,0+) 

3p 

aT-  =  ft  ui(s»°“) 
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so  that. 


Jim  W  u2(s.T(dT  +  i  t*U2(s,0+)  +  p  (s,0+)  +  p  (s,0) 

t* 

Jim  |  J  u2(s,T)dT  +  ^  t*U2(s,0-)  +  P?(s,0-)  +  p.(s,0) 

Z  -*“00  0  A  C 


P^s.O*)  =  P^s.O-) 


by  matching 


=  0  , 


and  by  Bernoulli 


U1 

P^(s)  +  —  =  constant 


U2(S,0+)  =  U2(s,0-)  =  U|(s) 

where  e  signifies  external  conditions. 

By  virtue  of  these  relations,  we  obtain 


(  l!  l2 

P2(s,°+)  -  P2(s,0-)  =  Jp2|  =  -jU/  «f(s,x)dT  -  /  u?(s,T)dx 

1  'wake  10  1  o  1 


-  <*!  - 


!  1 


166 

C347 lA/j  bs 


}  Rockwell  International 

Sctenca  Canter 

SC5055.21FR 


t^-*» 

[p2]  =  \  J  [uf(s,T)  =  U|(s,T)]dT 

=  I  f  fui  -  ue][ul  +  Ue]dT 

-co 

'll  ^ui  “  ue^uldT  +  TT  /  ^u2  Ue^dt  =  "  ~TT~  03) 

-00  -00 

where 

o  -  /  (Ug  -  =  momentum  thickness 

-oo 

Analogous  results  to  (13)  can  be  obtained  for  the  compressible  case 
from  (4b),  by  recognizing  that  the  significant  difference  will  be  folding  the 
density  factor  Inside  the  Integrands  of  s*  and  e  in  the  appropriate  definitions 
of  these  thicknesses,  and  by  employing  the  adiabatic  energy  invariant  rather  than 
Bernoulli  to  obtain  the  relationship  of  P1  and  Uj  in  the  outer  flow. 
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